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filtering. In this work the location of the eye of a 
hurricane is estimated from satellite observations. A 
comparison with a posteriori meteorologist's analysis 
was attempted. An adaptive gating scheme was employed 
in the filter to accommodate "manueuvers" of the storm. 

The observations were at random intervals and also 
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INTRODUCTION 



I . 

Satellite communications, especially digital, are well 
under way. One of the many observations that satellites 
provide are the meteorological. These give the location of 
a storm or a typhoon in terms of geological coordinates. It 
is possible to have images of the store or typhoon cloud 
cover . 

These observations occur quite randomly in time. 
Meteorologists like to present their forecasts equally 
spaced in time. 

The attempt, here, is to try to compare an adaptive 
Kalman filtering algorithm to estimate the location of the 
eye of a storm with optimum track values that are given from 
the meteorologist's analysis. The filtering process combines 
all the available measurement data with prior knowledge about 
the system. It produces an estimate of the location of the 
storm (latitude, longitude) in such a manner that the mean 
square error is minimized. The parameters, which are essential 
design elements of a Kalman filter, are the measurement 
noise covariance matrix, R, the excitation covariance, Q, 
the initialization covariance of error in the filter itself, 
P(l/0), and the transition matrix, 0. 

In most physical processes that one desired to track, 
many of these parameters change during the process of 
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tracking. The measurement noise associated with the 
observations can change if a different sensor is used, or 
a similar sensor obtaining measurements from some different 
geometry relative to the object being tracked. If the object 
being observed is acted upon by external forces, then the Q 
matrix should be changed to account for these external 
excitations. Most processes being tracked will change in their 
dynamic characteristics during the observation time and so the 
transition matrix ideally should also be changed. Further, 
the time between observations quite often occurs randomly 
in time. All of these things bring about the need to change 
the Kalman filter to adapt as the process changes. 

It is possible to change the parameters of the Kalman 
filter by sensing the error between the observation and the 
prediction from the track of what that observation should be. 

If a gate is established representing 95% of the normal 
random perturbations of the process, then when this error 
exceeds the magnitude of the gate, one can reasonably 
ascertain that the filter is no longer properly representing 
the observed process. In the work attacked here, real data 
was obtained from satellite observations and the qualitative 
observation error was established (PCN #). 

The error covariance matrix in terms of error ellipsoids 
along the track gives a measurement for the worthiness of 
the algorithm parameters. 
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II. KALMAN FILTERING TECHNIQUES 



In a linear, discrete system, the state and measurement 
equations are given by 

x(k+l) = 0(k)x(k) + fw(k) 



and 



^(k) = H(k)x(k) + v(k) 

where x is the state; transition matrix; _T is 

excitation noise matrix; H is the measurement matrix; w 
and V are the excitation and measurement noise correspondingly, 
assumed uncorrelated, zero mean white Gaussian: 

E[w(k) ♦w'^(j) ] = Q(k)6^j 

and 

E[v(k) •v'^(j) ] = R(k)6^j 

and 

E[w(k)] = 0, E[v(k)] = 0 

where Q(k) and R(k) are covariances of excitation and 
measurement noise. Now if x(k) is the estimated state 
value after the kth measurement and x(k|k-l) is the 
predicted value of the state before the kth measurement 
we have : 
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x(k|k-l) = f (x (k-1 1 k-1 ) , k-1 ) , where f is any function. 

The state error vector is defined to be 
x(k) = x(k) - x(k) 

and the predicted state error vector is defined to be 
x^(k|k-l) = x(k|k-l) - X (k) . 

The covariance of state error matrix is defined to be 
P(kjk) = E [x (k) -x"^ (k) ] , 

and the predicted covariance of state error is defined as 
P(kjk-l) = E[x(k|k-1) •x'^(k |k-l) ] . 

The state excitation matrix is defined by 

Q(k) = r(k)E[w(k) •w'^(k) ] -r^(k) , 

and the measurement noise covariance matrix is defined by 
R(k) = E [v (k ) • y"*" (k ) ] . 

The Kalman filter equations are: 

P(k+l|k) = 0(k)P(k|k)b'^(k) +Q(k) 

G(k) = P(k|k-l)H'^(k) [H(k)P(k|k-l)H'^(k)+R(k) ] 

P(k|k) = [I-G(k)H(k) ]P (k|k-l) 
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x(k|k-l) = _0 (k ) X (k-1 j k-1 ) 

^(k |k-l) = H(k)x (k| k-1) 

x(k|k) = x(klk-l) + G(k) [z (k) -z (k |k-l) ] 

The initial condition of P (error covariance matrix) and the 
Q and R matrices are the determining factors in the filter 
structure. For Q having main diagonal values greater than 
R means that we have greater uncertainty in the state 
estimate than in the observation. Thus the new state 
estimate is more dependent upon the new measurement and less 
related to prior estimates. The inverse is also true. For 
R having greater diagonal terms indicates that the new 
measurement are subjected to stronger corruptive noises, 
and so should be weighted less by the filter. The gains (G) 
are lower. The P is responsible for the initial transient 
performance of the filter. 
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III. ERROR ELLIPSOIDS 



The error covariance matrix in each stage of a Kalman 
filter process gives insight into the quality of the track 
occuring . 

The diagonal terms ^ 22 ^ variances of 

uncertainty in our knowledge of latitude and longitude. Their 
respective off diagonal terms are covariance between 
latitude and longitude. 

The square roots of the diagonal terms gives us the 
rms errors in out estimates of longitude and latitude. Having 
the definition of the structure we are dealing with (in our 
case the satellite observations) and its uncertainties 
(expressed by the PCN number-actually by the values of R) 
we can see how the K.F performs through its error covariance 
matrix. Expressing the P matrix in an ellipsoid of constant 
probability, one obtains a visual appreciation for the 
worthiness of the algorithm parameters. The representation 
requires that the errors are normally distributed. 

The joint probability density function is; 

^-l/2e'^(k| k-l)p"^ (k|k-l)e (k| k-1) 

where e(k|k-l) is the predicted state error vector. Setting 
the exponent equal to a constant value, we are going 
to have a curve which is an ellipse. This ellipse. 
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however, does not have its major and minor axis aligned with 
the coordinate system. Instead its axis (x', y'), comes 
from the following transformation: 

x' = XCOS0 + ysin0 

y' = ycos0 - xsin6 

where x, y are the latitude and longitude in our case with 
8 = 1/2 tan’h 

a - <T 
X y 

where cos(x,y) = P ^2 k-1) 

0-^ = P,, (k|k-l) and o-^ = P„„(k|k-1) 

X 11 I y 221 

The new variances are calculated by 
^2^2 

- a + a , . 

,,.2 ^ y ^ cov(xy) 

X 2 sin2 0 

- cr + cr , . 

,2 _ X y _ cov (xy ) 

y 2 sin2 3 

Incorporating the above equations, error ellipses are 
presented in subsequent figures with the satellite tracks. 
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IV. SATELLITE TRACKING-SIMULATION RESULTS 



Data from the Annual Tropical Cyclone Data for the 
Typhoon-Nelson appear in Tables 1 and 2 in terms of 
longitude and latitude. The best track data appears in 
six-hour intervals for twelve days and the satellite fixes 
(observed) in random time intervals. 

The data in the K.F algorithm parameters are: 



1 0 DT 0 I 

i 

0 10 DT 



0 0 1 0 
0 0 0 0 



r 



DT^/2 

0 

DT 



0 

DT^/2 

0 



0 



DT 



P(l/0) = 



10 ^ 0 
0 10 ^ 

0 0 
0 0 



0 0 
0 0 
10 ^ 0 
0 10 ^ 
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X(l/0) = 



0 I 

I 

0 

0 

0 



In Figure 1 we have a representation of the best track 
observed and the K.F track. An error ellipsoid at the 25th 
stage of the process appears. It seems that K.F follows 
the observed track more closely than the best track. This 
is due to the adaptive gating and Q relative to R. 

The gain history for G(l,l) fluctuates between 0.7-0. 8, 
never arriving at a stable value. The G(3,l) reaches a 
stable value of 0.12. 

The above appears in Figure 2 . The track for prediction 
appears in Figure 3 . 

The errors EBl and EOBl represent (YH-BLAT) and 
(LAT-BLAT) respectively. They, along with EB2 and EOB2 
appear in Table 3 and Figures 4 and 5 . 

Implementing Julian Time, we have a comparison in common 
time for 24 points only. It can be seen that in terms 
of the latitude error the K.F is close to the best track 
values . 

In an attempt for better performance the Q 
matrix was changed in the algorithm. Also the innovation 
errors, in terms of latitude and longitude, exceeding the 
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magnitude of the gate k-1) + R) resulted in changing 

the values of the gains. The above correction makes the 
filter more adaptive now and the errors EBl and EB2 appear 
smaller on an average. The representation of the above 
correction in terms of trajectories, gains, predictions and 
errors appear in Figures 6, 7, 8, 9, and 10 and Table 4. 

The computer program appears in Appendix A. 
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TABLE 1 



OBSERVED DATA (SATELLITE FIXES) 



TI ME 


LAT 


LONG 


IdO ^. 00 


3.70 


160. 90 


1 809. 00 


4.00 


160 . 10 


1813.00 


4 , 70 


157 . 70 


1900.00 


5.30 


155.50 


1903.48 


5 . 80 


154.30 


1906. 00 


5.80 


153.90 


1909.00 


6.00 


152.60 


l 9 t 6.33 


6 . 70 


150. 70 


1921.00 


6.80 


150. 00 


2000. 00 


7.20 


149. 60 


2003. 00 


7.20 


149.30 


2005. 1 8 


7.50 


149. 1 0 


201 2. 00 


7.00 


1 47. 40 


201 6.21 


7*20 


146 . 10 


201 8 . 00 


7.20 


146. 00 


2 02 1 .00 


7.20 


145 . 70 


2100.00 


7.30 


145.20 


2103. 00 


7.50 


145.00 


2105.06 


3.20 


1 44. 40 


2106. 00 


8.50 


144. 00 


211 2.00 


7.80 


1 43. 00 


2116.00 


8.00 


14 1. 90 


21 1 7. 51 


8.00 


14 1 . 70 


2200. 00 


8.30 


139.30 


2203. 00 


8 . 50 


1 39.00 


220 4 . 54 


8.50 


138. 10 


2206. 00 


b . 50 


138.00 


2212. 00 


8.60 


1 35.90 


2 ^ 1 6 . 00 


8.30 


1 34. 70 


221 7.40 


9.00 


134.20 


2300. 00 


8.80 


1 33.60 


2306. 00 


9.30 


132.30 


2312. 00 


9.40 


131.40 


231 7.28 


9 . 50 


130 . 10 


2321.00 


9.50 


129.60 


2400. 00 


9.50 


129.50 


2403. 00 


9 . o 0 


129.20 


2406. 13 


9.80 


128. 70 


2412.00 


9. 70 


128.20 


241 6. 00 


9.60 


127.30 


2418.00 


9 . 70 


127.60 


242 1.00 


9.90 


127. 40 


2500.00 


1 0.00 


127.30 


2503. 00 


10.10 


127. 10 


2506.01 


10.10 


126 . 80 


2509. 00 


10.30 


126. 70 


2512. 00 


10.30 


1 26.30 


251 8. 00 


10.40 


125.50 


2521 .00 


10.40 


125 - 1 0 


2600. 00 


10 .30 


124.50 
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TABLE 2 



BEST TRACK DATA 



time 


BL.AT 


BUONG 


1 306 * 00 


3.30 


160. 70 


1312.00 


4.50 


isa . 30 


lai a . 00 


4.90 


157 . 10 


1900. 00 


5.50 


155.50 


1906. 00 


5 . 90 


153. 90 


19 12.00 


6.40 


152.40 


1918. 00 


6 . 70 


151.00 


2000. 00 


6.90 


149. 90 


2006. 00 


6.90 


148. 70 


2012. 00 


7.00 


147 . 60 


201 a. 00 


7.30 


146. 50 


2100. 00 


7.50 


145. 50 


2106. 00 


7.80 


1 44.40 


21 1 2. 00 


7.90 


142. 90 


211 a. 00 


a . 10 


141.30 


2200. 00 


8.30 


139.50 


2206. 00 


8.60 


137. 70 


2212. 00 


8.30 


136. 10 


22 X 3.00 


9.00 


134.80 


2300- 00 


9.1 0 


133.50 


2306. 00 


9.20 


1 32.40 


2312. 00 


9.30 


131.40 


231 6.00 


9.50 


130.40 


2400 . 00 


9 . 70 


129.50 


2406. 00 


9.80 


128.80 


2412.00 


9.80 


1 20.20 


24 l 3 . 00 


9.90 


127.70 


2500. 00 


10.10 


1 27.20 


2506. 00 


10.30 


126.80 


2512.00 


10.50 


126.20 


251 a . 00 


10 . ‘♦O 


125.50 


2600 . 00 


10.30 


124.80 


2606. 00 


10.50 


124.00 


26 1 2. 00 


11.00 


123.50 


2613. 00 


1 1 .50 


1 22.^0 


2700. 00 


11.70 


122.20 


2706. 00 


11.80 


121.30 


2712.00 


1 1 .90 


120.40 


271 a . 00 


12.00 


119.50 


2800- 00 


12.40 


11 a , 50 


2606. 00 


12.90 


11 7,80 


2812.00 


13.30 


11 7.20 


2 aia. 00 


13 . 70 


11 6 .oO 


2900. 00 


14.10 


116. 10 


2906 - 00 


14.30 


11 5 . 30 


2912.00 


14.20 


115.40 


29 ia . 00 


14.00 


114. 90 


3000- 00 


13.90 


1 1 4. 50 


3006. 00 


13.80 


1 1 4.20 


3012- 00 


13.50 


1 14. 00 
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Figure 1 Trajectory of Storm Nelson 
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and velocity 




dHA 'HA 
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Figure 3 K.F Track and Prediction Ahead in Latitude 



TABLE 3 



ERROR VALUES IN TERMS OF LATITUDE AND 



.ITT .vi- 


r- ^ 1 


FIP 


FHl 1 


1P0^^4 P 


-O. Pp 


0 . r> 6 


-,i .PO 


7 'P 0 a. 


-O.Pl 


0.0^ 


-0 • PO 


7 POi-»^/i. 


- o . o t 


— O . T 6 


-0.10 


7 PO ;S7 P 


O . IH 


-0 .4^ 


o , 10 


7 P 0 r'i .1 a. 


O . PO 


-O .17 


0 • o 


7 P O K Q O 


-O . OP 


-0.6 7 


-o . 1 n 


7PO 60r> 


- O . P f't 


-0.11 


-O . PO 


7PO 7 1P 


o . 16 


-O .0 6 


0 . 70 


7P 0 70H 


-0 . 06 


0.01 


-0.10 


7 PO 7P O 


-O . 06 


- O . t C. 


0 .o 


7 PO 7P*^ 


-n. OP 


0 . P 6 


-0.10 


7PO 71P 


-O . P 1 


-0.1 ^ 


-o . PO 


7P0 744 


-o . op 


-0.1 1 


-o . in 


1 PO 7«1 O 


0.16 


-O. 1 1 


0.10 


7PO 7S/^ 


0.16 


-0.0^ 


0.10 


7P0 7r-. 1 


-O. 16 


- O . 0 rt 


-0 . PO 


7 P 0 7 7 4 


-o . OP 


-0.07 


O . 0 


7 P 0 710 


o . OP 


- O . 1 P 


-0.10 


7P0 71^ 


-O. 1 4 


-O .P<1 


-O . PO 


7 P O 7 O P 


-0.01 


-0.01 


- O . 1 ('I 


7 P O 7*0 A 


-0.00 


-0.07 


-o . PO 


7 P O H O u. 


-0.17, 


0.0 6 


-O .^0 


7P OP 1 O 


-O , 04 


. O 1 


r» . i‘i 
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— O . Or*. 
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0.0 



LONGITUDE 



0 . 60 
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-0.10 
-O . P(J 
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-o . in 
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-0. PO 
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n . o 
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Figure 6 Trajectory of Storm Nelson 



I 

i 



>1 

-P 



C *H 

o w 

•H o 







26 



Figure 7 Gains of Position and Velocity 
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Figure 8 K.F Track and Prediction Ahead in Latitude 



TABLE 4 



ERROR VALUES IN TERMS OF LATITUDE AND LONGITUDE 



J T T 


*^ri 1 


2 ^2 


2^ 1 


=■ j ^ ^ 


^20 2 


-0.21 


0 , ol 


.20 


0 T rv j 


7 J ^ O <4. .3 


- j -<,7 


0.^5 


- J .^0 


0 , :) 


y' ' L' 0 


‘ J . 1 


-0.12 


-0 . 1 0 


0 . 0 


72^ o72 


0.26 


- 0 . 7 


0 . 30 


-0 . 3 J 




0. i 3 


-0.21 


} . 3 


-0 , 


j :)0 


-0.12 


— 0 . 6 c5 


— • 10 


-7 . J 
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-0 . 2 3 


- 0 . ^ b 


— y . ^1. 7 


-0.3 -> 


770702 
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-0 , - 


0.7^ 


-0 , a,) 


7y^0 70.TI 
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-j . 
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7 20 7 2 j 


- 0 . 16 
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j . 1 0 
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— ) . 20 


— 0 . 2 J 


/ ^ 7^-4. 
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0 . 1 V 


- ; . 30 


0 . 10 


7^0 0 


0 . 1 0 


- 0 . ) 


C . 1 0 


-0.1^ 


7_i > 7 “n 


0 . 1 


-0 . 0 a 


) . 1 'I 


0 .0 


7^^' 7.;d 


-0. 14 


-0.71 


- j . 20 


0 . 0 


2 0 7 2 


-0. 00 


- 0 . 1 G 


:) . 0 


-0,10 


^^0 7--' 0 


-0, 04 


0.0a 


- . i 0 


ry 


■^20 73o 


— U . 2 0 


-V) . 1 I 


. 4. ,' 


-0 . 10 


720 7*^2 


-0,21 


- 0 , 0 7 


-0 . 1 0 


C . IJ 


"2'^ 7^^ *j 


-0.23 


-0 . C ^ 


- .^0 


0 , 0 


7 20 20 a 


-0 . lO 


- 0 . 6 


- . 2 7 


C • 1 u 


7 0 ci 1 0 
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Figure 10 Longitude Errors 



V. RANDOM TRACKING 



To examine the adaptability of the K.F, another storm was 
created. The equations for this simulated storm were: 

BLAT(S) = BLAT(s-l)-V -T 

X 

BLONG(s) = BLONG (s-1.) -V^-T 

LAT ( s ) = BLAT ( s ) + V ( s ) 

LONG(s) = BLONG(s) + V(s) 

where T=6hr, V =10°/24hr, V =5°/24hr and V = measurement 
X y 

noise (created by a random generator subroutine) . 

Implementing the above equations, a "true" and an 
"observed" track were created. These two data files were 
named BESTRACK and OBSERVED and appear in Table 5 and 6 . 

Simulating with the above new data, the K.F algorithm 
appeared to track well. Figure 11 shows the "true", 
"observed" and K.F tracks. 

Figure 12 indicates the gain history in terms of Gil 
and G13. This approach gave stable values of 0.27 and 
0.01 respectively after the 4th discrete point in time. 

After this time the gain does not vary any more. This 
means that the innovation error, ( z (k ) -x (kl k-1 ) ) is 
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weighted each time by the same quantity after the 4th 
observation. Having lower values in the diagronal terms of 
the Q matrix in comparison with R, in this case, means that 
we have greater uncertainty in the measurements (observed 
data) relative to the model uncertainties. So the gains are 
smaller and the filter no longer "tracks" the measurements 
closely . 

As far as the latitude and longitude errors are concerned 
it can be seen that EBl (YH-BLAT) and EB2 (XH-BLONG) are 
nearly zero and are smaller in comparison with EOBl (LAT-BLAT) , 
E0B2 (LONG-BLONG) . This shows the ability of the algorithm 
to follow the "true" values more than the measurement if the 
latter has been corrupted with noise. 

The error values and the plots appear in Table 7 and 
Figures 13 and 14 . 
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TABLE 5 



BESTRA.CK DATA 



TI ME 

IBOO . 00 

1806.00 

1812.00 

1818.00 

1900.00 

1906. 00 

19 12 .00 

1913. 00 

2000. 00 

2006-00 

2012.00 

2013.00 

2100.00 

2106. 00 

2112.00 

2113-00 

2200 .00 

2206-00 

2212.00 

2213.00 
2300-00 
2306-00 

2312.00 
2318-00 
2 ^ 00.00 
2406- 00 

24 12 . 00 
2413-00 
2500-00 
2506- 00 

2512.00 
2518-00 

2500.00 
2606- 00 
2612-00 
2618-00 
2700-00 

2706 .00 

2712.00 

2713.00 

2800.00 
2306-00 

2312.00 

2818.00 

2900-00 

2906. 00 
2912-00 

2913.00 

3000.00 
3006-00 



SLAT 

160.00 
157-50 
I S 3 .00 

152.50 

150.00 

147.50 
1 - 45.00 
1 42-50 
140-00 

137.50 
1 35-00 
132-50 

1 30 .00 

127.50 

125 .00 
122.30 

120.00 
X 17.50 
1 15 .00 
112-50 
1 10.00 

107.50 

105.00 

102.50 

100.00 

97.50 

95.00 
92-50 

90. 00 
87-50 

85.00 
82-50 

30.00 

77.50 

75.00 
72-50 

70.00 

67.50 

65.00 
62 .50 

60 .00 

57.50 

55.00 

52.50 

50 . 00 
4 7.50 

45.00 
42-50 

40 . 00 

37.50 



bi_ong 

3.30 

2.30 
1.00 
0.80 

- 0.20 
- 1 .20 
- 2.20 
- 3 . 20 
- 4 . 20 
- 5.20 
- 6.20 
- 7.20 
- 8.20 
- 9.20 
- 10.20 
- 11.20 
- 12.20 
- 13.20 
- 14.20 
- 15.20 
- 16.20 
- 17.20 
- 18.20 
- 19.20 
- 20.20 
- 21 .20 
- 22.20 
- 23.20 
- 24.20 
- 25.20 
- 26.20 
- 27.20 
- 23.20 
- 29 .20 
- 30-20 
- 31.20 
- 32.20 
- 33.20 
- 34.20 
- 35.20 
-36 .20 
- 37.20 
- 38.20 
- 39.20 
- 40.20 
- 41.20 
- 42.20 
- 43.20 
- 44.20 
- 45.20 
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TABLE 6 



OBSERVED DATA 



Ti'^e 

laoo, 00 
1806 » 00 
18 12 .00 

1818.00 

1900.00 

1906. 00 
19 12.00 

19 18. 00 

2000.00 

2006.00 

20 12 ,00 

20 13.00 

21 00. 00 

2106.00 

21 12.00 

2113.00 

2200.00 

2206.00 

2212.00 

2218.00 

2300 . 00 

2306.00 

2312.00 

2318.00 
2 * 100.00 
2 * 06. 00 
2 * 12.00 
2 * 18.00 

2500.00 

2506.00 
251 ^ i .00 

25 1 8. 00 

2600. 00 
2606- 00 

2612 .00 

2613.00 

2700. 00 

2706 . 00 

2712.00 

271 3. 00 

2800.00 

2306.00 
2812-00 

2818.00 

2900. 00 

2906.00 

2912.00 

2913. 00 

3000.00 

3006 . 00 



1_AT 



160.00 
157.30 
15*. 92 
152.71 
150.03 
1^7 • 6 3 

142.61 
139,99 
X 3T,26 
1 34 , 6 ;> 

I 32.39 

I 30 .3 1 

127,37 
X25 .19 
122.44 
120 . 19 

I I 7,40 
1X4 ,87 
X 12.92 
X 10,03 
107,65 
104 ,83 
102 .96 

99,90 
97,65 
95 .29 
92 .44 

90.5 1 
87,13 
84 ,34 
82 ,40 

80,00 
77,73 
7n . 0 8 
72 ,50 
69 ,96 
67,22 
64,79 
62,75 
59 ,83 

57.37 
54 . 72 
52,43 

50 .06 
4 7.91 

44.38 
42 ,53 
40 , I4 

37,2b 



LONG 

3,80 
2.60 
I .72 

1,0 1 
-0, I 7 
•1 ,07 
-2.25 
-3,09 
-4,21 
•5.44 
-6,58 
-7,31 
-7,89 
-9,33 
-10^0 I 
-11,26 
-12,0 I 
-13,30 
-14,33 
-14,78 
-16,17 
-17,05 
-13,37 
-18,74 
-20,30 
-21 ,05 
-21 .91 
-23,26 
-23,69 
-25,52 
-26,36 
-27.30 
-2B. 20 
-28,92 
-30, 12 
-31 ,20 
-32.24 
-33,48 
-34,41 
-34, 95 
-36,37 
-37.13 
-33, 48 
-39 ,27 
-40, 14 
-4 O • 79 
-42,32 
-43,12 
-44,06 
—45 , 44 
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Figure 11 Trajectory of Fictitious Storm 
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TABLE 7 



ERROR VALUES 


IN TERMS 


J I I r . iL " 


t h 1 


7 2062 ^ 


- 0.06 


7 20 6 ^0 


- 0.20 




- 0 . 12 


7 20 o4 2 


0 . 1 1 


720 o 4 d 


0.10 


7 ^ 06^4 


0 . 1^3 


720660 


0 • oo 


7 2 0 6 6 6 


0 . 1 1 


720672 


0 . 08 


72 0 67 6 


- U . 02 


7 2 0 4 


- 0.14 


72 ^^ 


-0. 1 5 


720o0o 


“ U . 02 


720 702 


-0 . o-'^ 


72070 a 


0. 02 


72071 ^ 


0. 00 


7 20 72 0 


0. 06 


7 20 72 D 


0 . 02 


7 20 77 2 


- 0.02 


■^ 207 Ja 


0.10 


720 74 4 


0 . 1 0 


720 760 


0.12 


7 20 70 6 


0 , 0 ^3 


720 762 


0.17 


720 76 6 


0 . 1 1 


7^0 774 


0 . 12 


720 7 --' 0 


0 . 1 v6 


720 7 Mo 


0.12 


72 0 70 2 


0.23 


720706 


0 . 10 


720 604 


0.02 


720610 


- 0.02 


720610 


- 0.07 


720622 


0 . 04 


7 2 0 62 6 


0 . 


7 20 6 3 4 


0.04 


720640 


0 . 0 1 


■"2 3 64 6 


- 0.07 


7 2 • ) 6 3 2 


- 0.12 


7 20 6-3 6 


— 0.04 


7 2 0 66 ^ 


- 0.08 


72 0^.70 


- 0 . 0^3 


72 0 67 


- 0.12 


72 3 6^.2 


- 0 . 1 1 


72 0 60 6 


- 0.08 


7 ^ 0 ^^ Q 


0. 05 


7 2 0 O 0 0 


0. 0 2 


72000o 


0.04 


7^001 2 


0.0 8 


7 2 0 V 1 8 


- 0 . 18 



LATITUDE AND LONGITUDE 





crO M 


FU82 


^0 


J . 0 


0 .0 


20 


— 0 .20 


- 0 .20 


1 J 


- ; .08 


-0 . 03 


1 0 


3.21 


0.2 1 


29 


0 .07 


0 .03 


1 J 


0.13 


0 . 1 J 


07 


-J . 05 


- 0 .06 


1 0 


0.11 


0 . 1 1 


0 7 


-0 . 0 1 


- 0 .0 1 


06 


- ) .24 


-0 . 24 


1 5 


. 38 


-0 . 38 


1 6 


-0 . 1 1 


-0.11 


0 3 


0 .3 1 


0.31 


Ot> 


-0.17 


-0.13 


0 1 


0.15 


0 . 1 0 


0 0 


-0 .05 


-0.03 


Olj 


0.19 


0.19 


0 ^ 


-0 . 1 0 


-0. 10 


0 2 


-0.13 


-0.13 


1 0 


0 .42 


0.42 


0 9 


0.03 


0.03 


1 2 


i> • 1 5 


0 . lb 


0 5 


-3.17 


-0.17 


1 7 


J . 46 


0 . 4 1. 


1 1 


-0.10 


, - 0 . 10 


1 2 


0.15 


0.13 


t 8 


> .85 


0 . ^0 


1 2 


-0 .05 


- 0 . Oo 


23 


J . b 1 


0 . 5i 


1 0 


-J . 3 2 


-0.32 


02 


-0 • 1 5 


- 0 . 1 0 


02 


-0.10 


- 0 . 10 


0 3 


L) . 0 


0 . 0 


04 


^.28 


0 . ?8 


0 5 


3 . 6 


0 . 0 0 


0 4 


0 . 0 


0 . 0 


0 1 


-0 . 0 4 


-0.04 


0 7 


. 2 8 


-0.23 


1 2 


-0 . ^ 1 


-0.21 


0 4 


J • 25 


0.23 


Ob 


— 0 .17 


-0.17 


0 5 


0 . 0 7 


0.0/ 


1 2 


-J .28 


-0 .2d 


1 1 


— w) .07 


-0.07 


0 6 


0.06 


0.05 


Ob 


0 . 6 1 


0.41 


02 


-0.12 


-0.12 


04 


0 . C 8 


0.08 


0 3 


0.14 


0 . 1 


1 8 


— .24 


-0.2^ 



OF 

-0 . 

-0 . 

-0 . 

0 • 

0 . 

0 * 

0 • 

0 . 

0 . 

-0 • 

-0 • 

-0 . 

-0 . 

- 0 . 

0 . 

-0 . 

o. 

0 . 

- 0 . 

0 , 

0 . 

o, 

0 . 

0 • 

0 • 

0 • 

0 . 

0 . 

0 . 

0 . 

0 - 

-0 . 

0 . 

0 . 

0 . 

0 • 

-C • 

-0 • 

- 0 . 

- 0 . 

-0 . 

- 0 . 

- 0 . 

-0 . 

0 . 

o . 

0 . 

0 . 

- 0 . 
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TABLE 7 



ERROR VALUES IN TERMS OF LATITUDE AND LONGITUDE 



J \ 1 ..It 


t 0 1 




r n . 1 


t‘u 0 2 


7 2 0 (.i ? ^ 


“U . OM 


- 0 O 


0 . 0 


0 . 0 


7 2 f * o 3 0 


- U . 2 0 


-0.20 


-J . 2 0 


-0,20 




- 0 . 12 


- 0 . 1 J 


-..00 


-0 . 0.3 


7 J ' ) n 4 2 


0. 1 1 


0.10 


3.21 


0.21 


^2 0 04 ^5 


0 . 1 0 


n . ^ ^ 


0.03 


0.02 


7^iOo*">h 


0 . 1 


0 . l J 


0 . 1 1 


0 . 1 1 


^ 2 > 0 


0. oo 


0.07 


- > . 0 5 


-o . f) ) 


7 2 0 f 1 o t > 


0. L 1 


P . ! 0 


0.11 


0 . 1 1 


720 0 72 


0 . OM 


0.0 7 


-J . 0 1 


-0 .O 1 


7jC o7o 


- U . (3 2 


- 0 . 0 J 


- ) .24 


- 0 . 2 


720 


-0.14 


- 0 . 1 * . 


. :ip 


-0 . 1 1 


7 2'^ J 


-0.1*^ 


- 0 . 1 o 


-u . 1 l 


-0.11 


720o‘)d 


-0 . 02 


- 0 . 0 i 


0 . o l 


0 . n 


7 2 0 7 0 2 


-U . 0*7 


— 0 . 0 t> 


— 0 .13 


-0.15 


720 70o 


U . 0? 


0 . p 1 


0 . [0 


0 . 1 y 


72071^ 


0.03 


-0 , 00 


-n . OP 


-0 . 0'^ 


7 2 » ) 7 > 0 


0 . 0 < > 


0 . 0 u 


;j . 1 o 


0 . l s> 


7 20 77l> 


0 . 02 


0.0^ 


-0 . l 0 


-0.10 


7 2 0 7^2 


-0.0^ 


-0. ')2 


-0.10 


-0.1 i 




0 . 1 0 


0.10 


J .42 


0.42 


72 0 74 4 


U . 1 0 


0 . 0 > 


C . 0 3 


0,03 


7 20 7S0 


0.12 


0.12 


u . 1 o 


0 . 1 S 


7 20 7':>o 


0.0 '*> 


0.05 


- } . 1 


-0.17 


7 2 0 7-2 


0.17 


0.17 


) . 4 f> 


n . 4 c 


720 7*» o 


0. 1 1 


0 . 1 1 


-0 . 1 0 


- 0 . 1 ij 


7, 0 7 7^ 


0 . 12 


0 . 1 2 


0 . 1 


0. 1 > 


7 2 0 7''. t) 


O . 1 M 


(3 . 1 'i 




) . ^') 


72 ) 7 Mo 


0.12 


0.12 


-0 .05 


- 0 . 0 o 


72 0 70 2 


U • 2 7 


0 . 7 J 


J . 5 1 


P . 5 1 


72') 7‘Ni 


0 . 1 0 


n . 1 0 


- J . ^ 3 


- ^ . 02 


7 jO M04 


0 . 0? 


0.02 


-0 . 1 5 


- 0 . 1 :> 


7 2 0 -J 1 0 


-0.02 


-0.02 


-a . 1 p 


- o . 1 . ) 


7 2 0 M 1 ( 'i 


-0.07 


- 0 . 0 J 


0 . 0 


f3 , 0 


720<-22 


0 . u 4 


v) .0 4 


^ • rj ‘ ‘ 


0 . ' » M 


7^0 02 M 


0 . 0^- 


0 . ) 5 


. L» 


0 . 0 c. 


7 20 0 < ^ 


0.04 


0 . .3 


i.' . 0 


0 . P 


7_ 0 j:* U 


0 . O 1 


0.01 


- . • .04 


-0.04 


3m 4 b 


-0.07 


- P . ^V 


- " . 0 


-0 . 2.-^ 


^ 2 0 ci ‘> 2 


- 0 . 12 


-P . 1 2 


- J . ^ 1 


-0.21 


7 2'^ o 3 M 


-0.04 


-0.04 


J . 25 


p . P > 


^2 3 rliM 4 


- 0 . OH 


-0.06 


-J . 1 


-0.17 


7 2 0^,70 


-0 . O'”) 


c 

* 

c 

1 


0 . 0 7 


0.0/ 


72 0^37 0 


- 0 . 12 


-0.12 


- ) . 2 ^3 


-P . 2 i 


/2 Oci *. 2 


-0. 1 1 


- 0 . 1 1 


-J . (J 7 


-0 .0 7 


7^ 0 


- 0 . 0 M 


-0.00 


O . 0 


0.0, 


7^0 


0 . 0 Tv 


p . O 5 


0 . ^ 1 


0 . 4 1 


720-000 


0.02 


0 . P 2 


-3 • 1 


-0.12 


7 ..'O ')<^ o 


0.04 


0 . 0 4 


0 .00 


0 . 0 O 


7<l 0 >12 


0.00 


P . 0 2 


0.14 


0 . 1 4 


7 20 V 1 H 


- 0 . 1 -M 


- 0 . 1 5 


- .' . ^ 4 


4 

'X 

. 

r 

I 
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72062 ^ 720664 720704 720744 720784 720824 720864 720904 

TIME 
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720624 720664 720704 720744 720784 720824 720864 720904 

TIME 



VI. CONCLUSIONS 



The K.F approach to estimate the storm's location appears 
to be very accurate. That comes from the comparison with 
the meteorologist's analysis results. Concerning also the 
fact that the latter was performed after the storm's occurence 
one can see the advantages of the K.F algorithm. During 
operation of the filter the actual residual sequence 

(k) (k| k-1 ) ] is compared to the gate [ Vp (kj k-1 ) +R (k) ] 
which actually is the square root of the residual covariance. 
Being a white Gaussian sequence, the residual is bounded by 
this gate. 

When the gate is exceeded the model is determined invalid 
within the filter and a modification in terms of Q and G 
takes place to adapt to the situation. At this point, and if 
the excess occurs in only one component of the vector residual 
process, one can further deduce that the measuring device 
generating the particular component is the source of 
difficulty (a sensor failure) . 

The error ellipsoids of the process also give insight 
into the filter performance in a more general case, referring 
to many sensors with a greater variety of uncertainties, the 
adaptive K.F algorithm could be a very advantageous approach. 
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APPENDIX A 
CCMPOTEfi AXGORIIHM 



Cl ‘iZ No ION HI (1 2, 12) , H ( 1 2 , 1 2) , F ( 1 2 , 1 2) ,G ( 12, 12) 

1 ,?IiIT (12, 12) , 

* Q(12,12), G31 (120) ,Pi<K (12, 12) ,?KKH1 (12, 12) 

1 , 3 1 1 ( IzO) ,G22 ( 120) 

DINENSICN IFEAD(10),IwHI'IE(10),Y(120),YH(120) 

1 ,XH (120) ,G4z (120) 

DIMENSION CEL(12,12),A(1z,12),a(12,12),D1(12,12) 

1 , 02 ( 12 , 12 ) 

DI'IENSIOM 2 (12) ,E (12) ,GE (12) ,XH? (120) , YHP (120) 

1 ,ZY (120) ,EX (120) 

DIMENSION DZLI ( 12, 1 2) , PHI ( 12 , 12) , II (1 20) , XKK'’'- (12) 



DIMENSION .'iAiiE (5) ,D (1^, 12) , X? (25) , Y? (25) ,KJDLH2 (120) , 
«h( 12,12) ,F (12,12),AI(12,12) ,ZKKi>!l(12),IK(120) 

1 ,?KK22 ( 120) ,P22 (120) , 

♦FT ( 1 2, 12) , ET (12,12) ,XKK ( 12) , XKXH 1 (1 ^) , ZXY (120) 

1 ,EOB1 (120) 

DIMENSION TIME (120) ,LAT (120) ,LONG (120) ,^11.1“ (120) 

1 ,SLAI (IzO) , 

♦ELCNG (120) ,BN1N 0 ( 1 20) , XK 1o (120) , YKZ6 ( 1 20) , PCN (120) 

1 ,£032 (120) 

IMIEGEF IXZ72 ( 100) , T K M48 ( 1 00 ) , I K X 2 4 ( 1 0 0 ) , 

♦ IKM 12 (100) , IK Mo (100) ,'IK?6 (100) , IKP12 ( 100) 

1 , XKP24 (IOC) , 

♦ TKP48 (100) ,IK?7z ( 100) 

OAT-A IQIT/1H£/,IY/'Y'//77;/’N’//IPH/' £'/ 

DEAL LA'I,ICNG 
I NT EG EH lYEMO /’Y’/ 
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I in ' G E H I Y E N 1 / ' Y ' / 

C 

C t-lZ= i'IC. OF C3SEEVELI V AL 'J EG ( S AT 3L L I T E ) 
C K2=NC. OF bZSl IHA^K VALUES 



;n=50 

dZ=50 

DC 132 1=1,4 
PKK(I, 1) = 1000. 

132 PKKiil (1,1) =1000. 

E ( 1 , 1) = . 00 00 0 1 

,M2,2) = . 00 0001 

C I'.EAC CE3ERVED VALUES 
C 

C HEAD (2, 1 1) (TI.'iE (I) ,LAT (I) , LONG (I) , PCN (I) , 1=1 , MZ) 

11 FORMAT (7X,?6. 2,F3. 1 , 1X,F4. 1 , 1X,f 1.0) 

FEAD (2, 1 1) (TIME (I) , lAT (I) , LONG (I) , FCN (I) , 1= 1 , NZ) 



C 

C READ dESI TRACK VALUES 
C 

READ (3, 1 4) (ETIME (J) , BLAT (J) , 3LONG (J) , EVIND (J) 

1 ,J=1,NZ) 

14 FORMAT (oX,F4.0,F4. 1, 1X,F4. 1,2X,F3.0) 

C READ (3, 11) {E1I3E (J) ,£IAT (J) ,BLOKG (J) ,?CN (J) , J=1,MZ) 

C 11 FCF.MAT (2X,4F10.2) 

C 11 FCR.1AT (oX,F4.0,F4. 1, 1X,F4. 1,2.(,F3.0) 



C 



ECHO VALUES * + * * 0 EMC V E NEXT THREE LINES TO ELIMINATE 
ECHO PRINT CHECK *♦#♦**«***:> 

WHITE (3,373) (TIME (I) ,LAI (I) ,LONG (I) , ETIME (I) 

1 ,3LAT (I) , 

* 3LONG (I) ,?CH (I) , 1=1, NZ) 

WHITE (3, 37 3) ( T Id E ( I) , L AT ( I) , LONG (I) , PCN (I) , 1= 1 ,NZ) 
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J73 



l'C£.:'iAT (2X ,7 F 1 0. 2) 

C ',.F I'JE (7, 37 j) (Ti:iE (I) , LA I (I) ,LO:iG (I) , E2I3E (I) 

1 ,HLAT (I) , 

C + DLCNG (I) , EL’INC (I) , 1=1, 'J2) 

C KaiTZ Of 37h) {II3" (I) ,LAT. (I) ,ICiiG (I) ,EIi;'jE (I) 

C 1 , SLAT (Ij , 

C * PCN(I), 1 = 1, N2) 

C 3 74 fCailAT (2X,6F10.2) 

C 



C THIS PaOGaAH CO.'IPUTES TFE FOLLOWING KALF.AN FILIIF 

C 

C - 1 

C G(K) = 2 (K/K- 1) *hT* (H + P (K/F- 1) 

C 



C ? (K/.<) = (I-G (K) *H) *F (K/K- 1) 



C 

C P(K/K-1) = CHI’!'? (K/K- 1) *PH11 -t-C 

C 

C 

L Q(I,J) DEFINES THE COVAFIANCE OF THE PEH SAHPLE FANDCH 
C EXCITAIICi'l Or THE FBCCES3 

C 



C 



C B(I,J) DEFINES THE 
C NHICB IS 

C 
C 

C H * ( I , J ) I S 

c 

c 

C II=K THE DISCRETE 



PAHOOH (GAOSSIAN) HEAS7EE.1ENT NOISE 
ADDED IC THE CdSERVAcLE SIGNALS 



THE IDEMIIT'i HATRIX 



POINT IN TINE, THE STAGE OF THE PRCCESS 
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c 



C i?KK{I,J) = i(K/K), (COV EFFOS AT K GIVEM K SA:''?LZS) 

C 

c 

C PKKK1(I,.J) = ? (K/K-1) , (CCV EEPOE AT K GlVZ.i K- 1 
CSAi''jELES) 

C N = NT.1B3H Of EOFJS, :i = N’U.TBE^ CF CCL.,C03. 

C , K C A l'( D d D A E 

C = NU.'-iBF.S OF I'lERATICNo C? THE FILTFF. 

C 

c 

NEVJ NEWS ^*’>'*****^4=^’!'^** + 

WRITE (6,7 17 1) 

7171 FCH.'IAI (/, 2 X /'****’*'******** H E W |g j WS** * ^***^ ^* *^*^* * * / , 



*2X, 


'THE STGR.1 


TRACK 


INPUT VALUES 


ARE A VIAL 


ABLE'/, 


=»2X, 


'THESE ARE 


ECHC 


FEIiNTEC IC THE 


TERMINAL 


AT THE'/, 


*2 X , 


'BEGINNING 


GF TH 


E PhCGSAM — A 


L 30 SEE L 


ISIING '/, 


^2X, 


*1? TOE Ei^CGRAO 


ENDS NORMALLY 


AN INPUT 


FILE 



1 :niL'/, 

*2X,'3E PRCDIJCED THAI HAY E£ PPINTEZ GZ USED 
1 F03 IN?UT'/» 

^zX,'t1SlNG THE SAHE FORMAT STATEMENTS 10 READ AS 
1 WESE'/z 

’>^2X,'USED TO WFITE Ctl UNIT 4. 33 *^ITH A FEN 

1 MOOIFICAIICNS Vz 

*2;(,'IHE BFANCH AFCHNn THE INPUT CAN BE USED. 

1 PRESENTLY'/z 

^2X,'A fllE CF THE INPUT DATA IS BEING PRODUCED AND 

1 WILL'/z 

*2X,'3E FOUND AS — K OUTPUT A-- Oi YOUR A - Cj.SK.') 
C 

C **** liiZ FIRST QUESTICN TC THE TER.MINAL — ASKS IF 
C AN INPUT FILE IS 10 EE USED 

C 
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1 



33 

C 

C 

C 

c 

C 

7777 

1234 

1235 
1016 
19 02 

10 

30 

40 



'.iBITZ (8,27) 
FC33AI (/2X , 



♦2X,'DG .N'CT USE ?A1 KEY 1C EXIT PROGBA'I UNLESS 
1 Y C U W A M I ' / » 

^2X'TO LOSE YCU2 INPUT FILE THAT IS PCINIED ON 



1 UNIT 4.'// 



V, 



*2X'EOR CHECK CUT ACTIVATE THE GO IG 883 STATEilENl 

1 CNV, 

*2X'1INE 135 TC END UITHCOT ERkOB'// 

»2X,'DO YOU AART TO DEVELOP THE INPUT FILE (Y/N)?') 



READ (2, 13) IANS 
FORMAT (A4) 

IF (IANS .NE. IYBN1) GO TO 9999 



FORMAT ( 1H 1 ) 

WRITE (8, 12 34) 

FORMAT ('ALL INPUT SHOULD BE FLOATING POINT EXi-EPT ') 
WBIIE(3,1z35) 

x-ORMAIC WHERE OTHERWISE SPECIFIED ',/) 

FORMAT (1H0,5X,5H X = , 13 ,/6 X, 5HGAI NS) 

WRITE (d, 10) 

FORMAT (/, 5X, • THE DISCRETE KALMAN FILIEF') 

W F 1 1 E ( 8 , 3 0 ) 

FORMAT (5X, ' ENTS? THE ORDER OE THE SYSTEM 

1 (UP TO 8) . ' ) 

R E A D ( 5 , 4 0 ) N 
FORMAT (I 1) 
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ND= 1 ^ 

D = 1 2 
LD=12 

DG 7d33 1=1, « 

Z (I) =0. 

CG 7393 J=1,N 
HI (I, J) =0. 

HI (I,I) = 1. 

?KI (I, J) =0. 

PHI (I, I) = 1 , 

K (I, J) =0. 

E (I, J) =0 
73^3 A (I, J) =0. 

« SITE {8 , '4 1 0) 

WRITE (8,7771) 

777 1 FGR3AT (5X, ' CO YO'J WAWT TC COMPUTE PHI S GAMMA 0 
1 LINE F?C;V , 

+ 1 X , ' A S 3 ? ' ) 

HEAD (5,7 77,;) TAN 

7772 FCPMATIAJ) 

IF(IAH .Ev.I?H)GC TC 7773 
IF (IAN . JE.IFH) GO TO 15 

7773 WRITS (8,7774) 

7774 FOHflAT (//,5X ,' ENTER THE A MATRIX') 

DO 7775 1= 1 , N 

DG 7775 J=1,N 
NFITS (8, 77 78) i , J 

77 76 FCKiMAT (5X, ' A (' ,11, ' , ' , 11 , ' ) = ') 

READ (5,7778) A (I, J) 

7778 tOSMAI (F10.3) 

7775 CODTINUS 
WRITE (9,7779) 
viSIIE (3,7779) 

77/9 FCPHAT (5X, ' THE A MATRIX •) 

CO 778J I=1,N 



46 



-h?V22 (-^,90) \A (I,J) , J=1 ,N) 

7760 V. t'lTE (6,00) (A (I,J) ,J=1,IJ) 

C >JF:iIS (6, 4422) 

C4-.22 FOnM AT (5X, ' EXIE5 TiiE CIME.^SION CF ') 

C 

C « IS CCI^iTEAlNED TC (1,1) HEPE o B SK 045000 0270 0 

C 

C 

C BEAD (5,442 j) 1 

04423 FCP'IAT (II) 

1=1 

«RIT£ (P,7731) 

7761 FC?.:iAT (5X, ' EMSS .HE 6 XATSIX ') 

DC 7732 1=1, N 

CO 7732 J=1,L 
WBirS (6,776 3) I, J 

778 3 ECPHAT (5X, ' E (' ,11, ' , ' ,11, ' ) = ') 

READ (5, 777e) £ (I, J) 

7762 CCNTIMOE 
.'('BITE (9, 7764) 

« BITE (8,77 64) 

7784 F0ETAT(5X,' THE E .'lATOIX ') 

DC 7785 1=1, N 

^.'SITE (9,90) (0 (1,0) ,0=1 ,1) 

7 78 5 h BITE (6, 70) ( E ( 1 , 0 ) , 0= 1 , L ) 

GG rC 9010 
15 X51TE (8,50) 

50 FO?OAT (5X, 'ENTEP IHE ELEOEHTS OF TEE IRAMSITON 
1 ;iArr<ix — pki') 

DO 1 1=1, N 
DC 1 0=1, N 
XBIIE (8,b0) 1,0 

bO FO'^P.AT (5X, 'PHI (',11,’, ',11,') = ') 

BEAD (5,70) PHI (1,0) 

70 FCFOiAT (F 10-0) 
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1 CCNriMUE 
i*eiIE(8,30) 

80 FOdriAT ( • 0 ' ,5X, • InE Oril X.VTHIX {T.XASSII ION XAIIIX)' 
DC 2 1=1, N 

2 SRiTE(b,i-0) (PHI (I, J) , J=1 ,N) 

90 FOfiaAT { 1P7E1 1.4) 

91 format (1P7I1 1. 4) 

C 

C GO ro 883 

1000 WFII E (3,100) 

100 rCRXAT (5'C, ' CO YOU WANT TO CHANGE ANY ELE'IENT OF IH 
1 NAIHIa?') 

900 HEAD (5, 1 10) IAN 
no FCFNAT(A3) 

If (lAN.Zg. IZ) GOTO 19 
IF (IAN- NS. lY) GOTO 1000 
-PITZ (3, UO) 

120 iOENAT (bX, • IvHICH ELSZENT OF THE MAIFIX 00 YOU NANI 
1 TG CHANGE?',// 

+ 5X,'SN"iEH AS IJ; WHERE I IS THE RCii AND J IS IHE 
1 COIUNN.') 

BEAD (5, 130) 1,0 
FOROAT (211) 

WRITE (8,50)1,0 
READ (5,70) PEI (1,0) 

WRITE (8,30) 

DO 3 1=1, N 

wRIIS (3,90) (PHI (1,0) ,0=1 ,N) 

WRITE (8,140) 

FCHil AI (5X, ' AN Y CTKER CHANGES?') 

GOTO 900 
WRITE (9,14 1) 

F020Ar(’1',5X, 'DISC BETZ TIME', 



130 

3 

140 

19 
14 1 
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+ IX, • ,<AL':\N rilTZ^ PPCGPAX') 

(9, li^2) (:, A/.E (Ij ,1=1,5) 

142 FCP/.AT (&X, ' FrCbLEi-i I D £ N 1 1 r IC ATIC , 5X , D A 4 ) 

'a = II£ (9, 14 3) 

1 43 FG&.IAT { ' 0 ' ,73 { ' * ') ) 

WHITE (9 ,30) 

DC 1111 i= 1 , N 

1111 WRITS (9, 91) (9HI (I, J) , J= 1 ,N) 

WRITE (3, 4 1C) 

WRITE (3, 150) 

150 FOPd AT {5X, ' ZWTEP. IHE DIMENSION CF FAMDOd INPUT 
1 VECTOF (W) ') 

READ (5,40) I 
WHITE (8, IbO) 

loO FORMAT (5X, 'INTER THE EIEXENIJ OF TUE D IS TFa C 0 I IC N ' , 
+ IX, ' dATF.IX — OAfl'dA. ' ) 

DC 4 1=1, N 
DO 4 J=1,I 
WRITE (8, 17 0) I , J 

170 FORdAl (5 X, 'GAdEA ( ' ,I 1 , ' , ' , I 1 , • ) = • ) 

READ (5,70) CEL (I,J) 

4 CONTINUE 

W R I T E ( 8 , 1 8 C ) 

loO FCRXAT {' O' ,5 a, 'THE GAi'ldA MATRIX (D 1ST R 1 3 UI ION 
1 d A i R I X ) ' ) 

DC 5 1=1, I'l 

5 WRITE (8,90) ( C EL (I , J ) , J = 1 , L ) 

1100 wRirE(a,iOG) 

9 10 READ (5,110) IAN 

xF (IAN . SQ. IZ) GOTO 29 
IF (IAN. NS. lY) GOTO 1 100 
WHITE (S, 12 0) 

READ (5, 130) I , J 
WRITE (8, 17C) I, J 
?EAD(5,70)CEL(I,J) 



49 



6 



rt jITi] (8, 16 0) 

DC b I=1,N 

;.'?ITE (8, 90) (DEL 1 1, J) , J = 1 ,L) 
rtSITE (3, 1UC) 

GOTO 9 10 

2 9 KBITS (9, 130) 

DO 1112 1= 1 , .N 

1112 K3ITi:(9,91) (CEL(I,J) ,J=1,L) 

WPIIE (8,4 10) 

90 10 K«ITE(3,300) 

3J0 FOR:1»\T (S/I, 'ENTER THE ELEMENTS OF THE CCV OF 'i • ) 
C 

C NOTE K IS CCNSTIifllirEO 10 (1,1) s) KK 17 10 & 4450 

K { 1 , 1 ) = . 0 0 0 0 0 1 
C 

7c BO 21 1=1,1 
CO 21 J=1,L 
KSTIZ (8, JIG) I, J 

310 FORM AT (5X , ' , I 1 , ' , • , 1 1 , ' ) = • ) 

READ (5,70) rt (I, J) 

21 CONTINUE 

77 WHITE (8, 32 1) 

WRITE (9, 32 1) 

321 FORMAT ('O' ,5X, 'THE CCV OF W) 

73 DO 22 1=1,1 

22 WRITE (3, 90) (S (I, J) ,0=1 , L) 

1400 WRIT£(8,100) 

940 READ (5, 1 1 0) IAN 

If (IAN. E;. IZ) goto 69 
IF (IAN. NE . TY) GOTO 1400 
WRITE (8, 120) 

READ (5, 130) 1,0 
WRITE (3, 310) 

READ (5,70) W (I , J) 

N R 1 1 E ( 3 , 3 2 1 ) 
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i^C 

2 3 W P 1 1 E ( d , 9 0 ) { ;; ( T , 3) ,3=1,1) 

VeiTE (8, 14C) 

Goto 9-40 

89 DC 1115 1=1,1 
1115 'aRITZ (9,9 1) (W (1,3) ,3=1,1) 

WRITE (d, 190) 

190 FCRDAT (5X, 'ENTZP THE CFDEE OF H, I.E. 3=?') 

READ (5,40) M 
WFIIE (8, 195) 

195 FOPZAT (5X, ' ENIE'1 IRE EL£FiE3TS OF THE CdGFRVAIICN 
1 .'1AT3IX--5. ' ) 

DC 7 1=1,:-! 

DC 7 3=l,il 
WHITE (3,20 0) 1 , 3 

20 0 5CH3AT (5X, ' H ( ' , 1 1, ' , ’ , I 1 , ' ) =' ) 

HEAD (5, 70) H (1 ,3) 

7 CONTINUE 

WRITE (8,2 10) 

210 FORMAT (’O' ,5X, 'THE H OATSiX (0 jS EP. V A T 1C N MATRIX) ') 
DC B 1=1,3 

3 WHITE (3,90) (H (1,3) ,3=1, N) 

1200 WHITE (8,100) 

920 HEAD (5 , 1 10) I AN 

IF (IAN . EQ. IZ) GOTO 39 
IF (IAN. NE. lY) GOTO 1200 
WRITE (8, 120) 

HEAD (5, 130) 1,3 
WRITE (8,200) 1,3 
READ (5,70) B (1,3) 

WRITE (3,210) 

DO 9 1=1,3 

9 WHITE (8,90) (H (I,J) ,3=1 ,N) 

WRITF (8, 140) 
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GCiO 920 
J9 0 9 ITH (9 , z 1 0) 

DC 111J 1=1, 

1113 '..5112(9,91) ( H ( I , .T) , J = 1 , N ) 

WRIT2 {'5, 193) 

WRITE (3,270) 

2 70 FOR.'IAT (5X, 'ENIEP THE ELE.*1E.NTS OF IHE 3EA30RE iESI 
1 NOISE’, 

+ 13 , 'COVARIANCE MAISIE — P ') 

7 2 DC 116 1=1,. 'I 
DC 116 J=1,,M 
WRITE (d, 280) 1,0 

28 0 ?0.-v;'!AT (5X, 'rv {' ,11, ' , ' ,^1, ' ) = ') 

HEAD (5,70) E (1,0) 

116 CCNTIMOE 

WRITS (8, 2'JC) 

WRITE (9,2JC) 

290 FCEMAT (' 0 ' ,5X , • IHS R dATTIX {idEASiJES.1 E NT • , 

+ 1X,'MCISE CCVAPIAOCE EAIRIX)') 

74 DC 17 1=1, H 

17 WHITS (3,90) (P (1,0) ,0=1,.1) 

1 300 «FITE (3, 100) 

930 READ (5 ,110) IAN 

IF (IAN. SC. 12) GOTC 79 
IF (lAN.NE. lY) GOIC 1300 
WRIIE (8,120) 

HEAD (5, 130) 1,0 
'WRITE (8, 230) 1,0 
HEAD (5, 70) x5 ( I, 0) 

WRITS (3,290) 

DC 13 1=1, N 

13 WRirE(3,90) (P (1,0) ,0=1,E) 

..RITE (3, 140) 

GOTO 930 

79 DC 1114 1=1,3 
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1 1 ]'4 



UKIIE (9,9 1) (? (I,J) , J=1 ,ii) 

wraiz (3, 4 1 G) 

(3, 3jC) 

:;30 FC5;iAT {//,Sa , ' r h K 'i 1 ( 1/0 ) 

42 DO J4 1=1, ti 
DC 34 J=1,N 
«RIT3 (8, 34C) I, J 

340 FORMAT (5X, 'PFK3 1 { ’ , I 1 , ' , ' , II , ' ) = ' ) 
RZAD (5,70) rKKcl 1 (1,0) 

34 CONTIWUF 



') 



•*3 WHITE (3, 351) 

WHITE (9 , j5 1) 

331 FC33AT (' 0 ' ,5X, ' ?HK3 1 ( 1/0) ') 

44 DC 35 1= 1 , N 

3 5 ii HITE (0,9 0) (FKKX1 (I,J) , 0=1 ,M) 

1500 0HiTE(3,100) 

950 HEAD(5, 110)1 AN 

IF (lAN.E^. IZ) COTO 51 
IF (IAN. :iE. II) GOTO 1500 
WEIIE (8, 120) 

?E AD (5, 1 30) 1 ,0 
WHITE (3, 340) 1,0 
HEAD (5,70) ?KK31 (1,0) 

WFITE (8,35 1) 

CO 3d 1=1 , N 

3b •„?. lTE (d, 90) (OKKOI (I ,0) ,0=1 , J) 

WFI :i (8,140) 

GOTO 950 

5 1 DC 1116 T= 1 , N 

1116 WHITE (9,91) (FKK.n (I,J) ,0 = 1 ,4) 

WHITE (9, 14 3) 

WHITE (8,246) 

240 fCP:‘!AT (5X, ' ENIEF, THE N0.1EEH Of THE 9CINTS TC 3E 



1 FEHr'CHMED. ' , 

+/,6X, ' (<10J) 7HIS Id AN INTEGER VALOE') 
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“E AD (5 , 247) 

2 47 rC;MA:(I2) 

IT= { 1) - 14J0) * (.0 1) 

IS= (2) -140 J) * (.01) 

1=24*11 

S=24*IS 

T=TId£ (1)-IT*100-14Q0 + I 
S = TI.1E i2)-IS*1D0-l400 + S 
DI=S-T 

27d3 continue 
Jh= 1 

4? II£ (d, 78 65) JK 
7365 rCa:-] AT (14) 

WHITE (8,51 1) ND , i'll , ID , HN , DI 

51 1 FOfeMiiT (2X, 2fiN=, 15, 5X,2H:1 = , 15, 5X, 2HL=, I5,5>:, duND 
1 ,15 ,5X, 3HMD, 15, 

*5X, 3tiLD=, 

8I5,5X,3HNN=,I5,5X,3KDT=,F10.4) 

WHITE (3, 533) 

533 FORMAT (/' MATRIX R ') 

C 

DO 3017 1=1, M 

3017 «?ITF(8,90) { K { I , J) , J= 1 , M) 

WHITE (8,544) 

544 FORMAT (/' MATRIX Q •) 

DC 3018 I=1,N 

30 13 5,’HII£ (8 ,90) ( Q ( I , .1) , J= 1 , » ) 

TRITE (3,555) 

555 FORMAT (/• MATRIX HKHMI') 

DO 3019 1=1, N 

3019 WRITE (3, 90) (PXKM 1 (I , J) , J= 1 , N) 

C IF (IANS. NF . I?U) GO TO 6 73 9 



54 



1-1= (ZlAl ( 1)-1400) 'I- (.Tl) 
15= (li;i£ (2)- IlOO) * (.01) 
'1 = 24*11 



S= 2‘4'i'IS 

I = TI'1S (1) - IT* 1 JO- 1 400 + 1 
S=Tl;iE (2) -13*100-1400 + 5 
5'I = 3-T 



JK = 2 

(8,7365) JK 

9753 CCNIIMUE 

y5ITt{5,7a76)CT 

7 878 fCR -lAT (5X, E 10, 3) 

CALL PHIDEL (DT,:i,L,A,5,PriI, DEL, 0 1,02,NO,.iO,LO) 

6789 iJRIIE (d, 666) 

6o6 FOFMAT(/' PHI ') 

DO 30^0 1=1, N 

30 2 0 WeilE (8, 9 0) { PHI (I , J) , 0= 1 , N) 

nP.ITE (d, 7 7 7) 

777 FOR3AI (/ ' DEI ' ) 

DO 3021I=1,N 

302 1 'JP.ITE (8,90) ( DEL (1 , 1 ) , 0= 1 , L) 

CALL TRANS (CEL,N, 1,CELI,HD,:iD) 

CALL PRCC (DEL, DELI , M , 1 , M , Q ,NC,:’IC,LD) 

CALL CO MSI (X (1, 1) ,2,N,N,Q,ND,.ID) 

LKITE (3, 544) 

DO 3025 1=1, N 

302 5 WRITE (8,90) (Q (I , J) , J= 1 , N ) 

4RITE (8,444) 

444 FOR!1AT(/' H ') 

DO j026 1= 1,3 

302o W?.irS(8,90) (H (1,0) ,J = 1/M) 

J K= 3 

WRITE (8, 7865) JK 
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rt piTZ 3 1 11) , -i, L, .n, '1C , L D, :;n, Di 

5111 ? C P. ■I A I ( 7 1 4 , Z 1 C . 4 ) 

DO 3 21 1=1, N 

3327 St:iTE(4,«0) ( A (I , J) , J= 1 , N) 

DO 5328 1= 1 , i\ 

3 328 «PirF(4,90) E(I/1) 

DO 5329 1=1, K 

5329 a5ITE(4,90) ('2(1, J) ,J=1,N) 

DO 533C 1=1, K 

5330 i:SITZ(4,90) { E K Kd 1 ( I , J ) , J= 1 , N ) 

DC 6327 1=1, y 

b3z7 .;EITF(4,9 j) ( a ( I, J) ,.J=1 , Li) 

DO 6927 1= 1, 

6927 i;siTE(4,90) ( n I (I , J) , J= 1 , N ) 

DC 6328 1=1, d 

6328 4[iITE(4,90) ( E ( I , J ) , J= 1 , .1) 

URITZ (3 ,7777) 

99 99 CGliTIliOE 

IF (IANS. EC. I SEN 1) GO TO 7 23 4 
IaIRITE (8, 1928) IAN’S 

1923 ~C?dAT(5X,' XX x XXX .TN XXXX X X *,?'4) 
JX=35 

ViFIIE (3, 7865) JE 

5EAD (4, 51 1 1) -N , d,L, ND ,dC,LD , ;iN, D I 
i] K = ) 3 ) 

Vi BITE (c, 7665) JK 

wBITE (8, 51 1 1 ) N , d, I, ME , :1C , LD, N t\. Cl 
JK = 4 

WRITE (3, 73c5) JK 
DC 723j 1=1, N 

7235 READ(4,91) (A ( I f 3) , 1 , ^i) 

JK = 432 1 

WRITE (8,7865) JK 
DO 7236 1= 1, N 
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12 jc 3ErtD('-f,9 1) E(I,1) 

DC 7^38 i— 

7z3d FEAD(4,'31) ( Q ( I , J) , J= 1 , N ) 

DO 7239 1=1,. N 

FEAD (4,91) (FKKill (I,J) , J=1 ,N) 

7239 W3ITZ{8,90) (? KK3 1 (I , J ) , J= 1 , N) 

DC 7328 1=1, a 

7326 P2AD{4,91) (U (I , J) , J= 1 , N) 

CO 7329 1=1, N 

7329 FEAD(4,4i) ( HI ( I , J ) , J= 1 , N ) 

CO 732 7 1= 1, M 

7327 3EAD{4,91) ( 9 ( I , J ) , J= 1 , H ) 

IT= (TIrtE (1) - 1500) * (.C 1) 

13= (113’=- (2) - 1500) * (.0 1) 

1=24*11 

S=24*IS 

I = IT3S (1) - II’*' 10 0-1500-t-I 
S = TI3S (2)-IS*10 0- 1500 + S 
DT=S-X 

9253 CCN'IIMUE 

X (1, 1) =.90 00 01 
JK = d 

W FITE (8, 73 6 5) JK 

CALL TAAiiS (CEL , 9, 1 , DELI , ND , 30 

CALL PROD (DEL, CELT, 9, 1,9,9 , HE , HD , LD) 

CALL CCN3T (a ( 1 , 1) , 9 , 0 , N , 9 , MO , 3D ) 

CALL PHIDEL ( DT , 'I , L , A , B , F HI , DE L , C 1 , D 2 , A C , 9 D , I D) 
7234 CCNriMUH 
DDT=DT 

a'EITE (3,91Ez) 

9182 fCa:-lAr(5X, ‘ ICO ACE AT DO 2222 ' ) 

C DLDLDLDLLDCLDOLDLCLDLDLDCLCILCLDLDLCLDIEDLD 
EB= 1500 + 11* 100 
KKE= 1 



57 



UO 1222 K = 1 , 

c 

C CALLS SUDHOJII.'. E TC CALCULAiE JULIAN IINE 

C ?c'^ tVEFY SiG?..'l POSIIICN AND EVERY HOI.ISS 

i'ii:iE=iNi ( :i:iE (K) ) 

IDAY = I'ri.lE/100 

IHCUh = ITi:iE- (IDAY^IOO) 

IF (IHOUR.EC. 0) IHC0R=24 
aCD6 = ;iOD (INCUR ,6) 

CALL JUTlaE (ITI>1E,JULHF) 

C IE {i10D6. EQ-0) GOTO 198b 

;iTia£=i:ii (ee) 

MDAY = '1TI NE/IOO 

:-lDAY = :iDAY* 100 

rfHOUH=aiiaE-HCAY 

IF (U HOUR. SI. 24) ariW£=IIia2 

CALI JUIIME (aTiaE,aJULHR) 

193b CCNTINUE 

'/iRIIE (8, 19 8S) JULnR,IIiaE 

1I8S FOR'IAT (///, ‘ JULIAS HOUR IS ’,19,' , ACTUAL Tia£ IS; 

1 ’,16) 

UEilTH (8, 19S4) aiias,ajuiHE 

1984 FOiiaAXC aCCULA 6 Ti:iE= ’,15,’ , COR 5 E 3PONDI NG 

1 JULIAN LE= ' ,19) 

1985 CGNTINOS 

C *** END JULIAN TiaE nCUIINE 

C 

C CAICUI. AIE aODiJLA 6 JULIAN TIJS f C F rKP,:KJ,XK? C 

, xsa ; 

r 

C *** .'lOCULA 0 EOn JULIAN IIME 

C IE (MHCUF.SC. 18) GOTC 3187 

C 

IKN72 (K) ^dJULUE-72 
IK1148 (K) =d JULHFs-43 
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'IK.^2 4 (ft) =2IjniHI^-24 
IK.’IU (K) 

'ih'/iG (K) =.ijclh:;-c 

IK (K) =1J'JLnt+6 
TKF12 (K) =.1.]UIKK+ 12 
iy?24 (K) =3JULHH+2^ 

IKPag (K) =« JULHF + 48 
IK?72 (K) =.iJUIf:B + 72 
3 1S7 CCN'IINU': 

C END jrjIIAN ri-IE POIJTINE 

C 

12 (PCiMf'") -5. NH . 0) GO TO 3133 
5 {1,1)=-25 
H (2,2) =.25 
GO TO 3134 

3133 IF (?CN (K) -3. NE. 0) GO TO 3135 
5 (1, 1) =.0c25 

R (2, 2) =.0625 
GC TO 3134 

3135 IF (?CN (K) -2. KE. 0) GC TO 3134 
B ( 1 , 1 ) = . 0 j 1 z 
5 (2, 2) =. 03 12 

3134 CG3TIMUS 

4EITE (8, 5445) lI.-iEIK) 
d 44 5 FORMAT {//////3 OX , 5 H T I M E= , r 1 0 . 4) 

■sRITE (4,5445) IIE5 (K) 

WIUIE(3, .313) K,3£,B (1,1) ,C (1, 1) 

9313 FORMAT (3X, 'K=' ,13, 5X, '85 = ' ,?B. 2,5 X, *£=' ,?7.-.,5:(, 

1 ' Q ( 1 , 1) = 'F 10 . 4) 

XPITE (4,93 13) K,B3, 5 (1,1) (1, 1) 

« RUE (8, 3 1 3) pen (K) ,DT , X ( 1 / 1) 

3 13 FORM AT (3 X, ' n ( K) = ' , F 1 0 . 4 , 3 X , ' ET= ' , F fc . 2 , 1 0 X , ' X ( 1 , 1 ) 

1 =',F10.4) 

WRITE (4,3 1 3) PCM (K) , DT , 'w (1,1) 
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DO 3i2'j 1=1,:'; 

DC 3129 J=1,H 
3129 Q (I, J) =0. 

Q{1,1) = {DC=^*4/^) (1, 1) 

Q(2,2) = (ET=^’;‘4/4) (2,2) 

Q(3,3) = (2 ,2) 

Q (4, 4) = DI**2*\': (2,2) 

/JRITZ (3,79 9) 

7 99 F0D:1AT (/' :'1AIRIX Q ') 

DC 3123 1=1, N 

3123 «SrTE (3,90) (C ( 1 ,0) , J= 1 , N) 
w (1 , 1) =.00 0001 
(2,2) =. 000001 

Call PHIDEL (E1 , i'l, I. , a , E,PHI , del. Cl ,D2, ad , .'ID , LD) 
93I1E (3,979) 

979 FG3:'1AI(/» Pril ') 

wailE (4 , 979) 

DC 3 97 9 1=1, li 

ARITE (4,93) (?.3I (I,J) ,J=1 ,N) 

3 5 79 KEITE (8,90) ( E h I ( I , J ) , J= 1 , N ) 

CALL GAi;> (PJ<K,?KK:1 1, ;,S, ? 31 , H , 4 , , 0 , H i , w D , 3 D 

1 /LD,K) 

WRITE (4, b5o) 

WRITE (3,656) 

6 5d format (/' PKK ') 

DC 3023 1=1, N 

WRITE (1, 90) (P:'<K (i,j) ,J = 1 ,:i) 

3023 ;iEITE{4,90) (PKK (1,0) ,0=1 ,:i) 

CALL FRCD (? I; I , X K K , N , N , 1 , X.\ K M 1 , ND , .X D , LD) 

CALL FECD (H, XKK;'! 1 , N , :i , 1 , ZKE.II , ND, XC, :1D) 

w^^IrE (3, 33 10) 

38 10 FORMAT (/' Z.^KMI ') 

WPIIE (3,90) (ZEKM 1 (J) ,0=1 ,0) 

WRITE (4,90) (ZKK8 1 (0) , 0= 1 , M) 

WRITE (d, 36 19) 
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dS 19 



FO ? :•! aF (/ ' l.A T ( , LG .13 ( K ) ' ) 

C W2I1H (3, 90) I,\T (K) ,LC:i3 (A) 

Z (1) =IA'r (K) 

L (2) =LO!.G (‘^) 

'..’HUE (4,08 11) 

'.^RITS (4, 90) (Z (J) ,J=1,i1) 

'.i?ITE (3,38 11) 

WRITE (3,90) ( Z (0) , J- 1 , *1) 

CALL SUB (2,ZKa'I 1 , 3 , 1 , E , N D , 3D) 

WRITE (3, =445) TIMECK) 
c 'nPiTE (3, ea 10) 

WHITE (8, 3029) 

3029 E''ORMAT(/ ' E IJ 

WRITE (3, SO) { S (J) ,0=1 ,11) 

If (K. LE. 1 ) GO TO 2204 
GAIE= (?KHi1 1 (2, 2) +Z ( 1 , 1) ) **.5 
I? (ABS (E (2) ) -GATE .11.0.) GC 10 2203 
G (2, 2) =0. 5+ ( 1. 2+G (2, 2) ) 

W (2,2) =1O0C0. ♦ 'W (2 , 2) 

G (4, 2) =0. 5=!“ (0. 33 3 + G (4,2)) 

C FK.CI 1 ( 1, 1) =2*FKcCOl (1 , 1) 

C PKKM1 (2,2) =2*F'KKZ1 (2,2) 

C RKKitl (3, 3) =2*PKKi1 1 (3, 3) 

c PKKFi 1 (4,4) =2*PKK.n 1 (4 , 4) 

C WRITE (4,9 192) GATE, E (2) 

WRITE (8,9 192) GATE, E (2) 

2203 CCNTIWUE 

GATE = (PKKd 1 ( 1 , 1) +?. { 1 , 1) ) * + . 5 

IF (ABS (Z (1) ) -GAT E .11.0.) GC IC 2204 
G ( 1 , 1) =0. (1 . 2+G (1 , 1) ) 

'W ( 1, 1) = 10000. * « (1 , 1) 

G (3, 1) =0.5’> (0. 333+G (3,1)) 

C PKKdl (1,1) =2*FKKM 1(1,1) 

C PKK11 (2, 2) =2*?KKil1 (2,2) 
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c 



c 



c 



yKK.11 (3, 3) =2*PKK.11 (J,3) 
FKK.1 1(4,4) =2*e.KK,1 1 (-.,4) 
W ?IT E ( 4 , 9 1 9 1 ) o AT Z , 2 ( 1 ) 
(J,9 19 1) GATZ,E (1) 



9 191 rCr 1AI (9X, ' EFFCT JT GATE. GATZ= ' , F 1 0 . 4 , 9 K , ' Z ( 1 ) 



1 ,E 10. 4, 'XXX') 




1 , F 1 0 . 4 , ' X X X ' ) 



C 



c 



c 



c 



c 



G1 1 (K) =G { 1 , 1) 



C 



C 



c 



931 (K) =G (3, 1) 

DC 3022 1= 1,11 

k’PITZ (8, 90) (FKKM1 (I,J) , J=1 ,N) 



30 22 KRI?E(4,90) ( ?K KM 1 ( I , J) , 0= 1 , !l) 

2294 WRITE (4, S9) 

WRITE (3,99) 

9 9 FORMAT (/' M. IIP IX 0 ') 

DO 3024 1=1, M 

WRITE (3,90) (G (I, J) ,J=1 ,M) 

30 2 4 WRITE (4,9 0) { G ( I , J ) , J= 1 , M ) 

WRITE (S,90) (ZKKM 1 (J) ,0=1 ,M) 

WRITE (3,38 11) 

WRITE (4,38 1 1) 

83 11 FORMAI(/' Z ') 

rtRITE (3, 90) [Z (J) , J= 1 ,[']) 
rtRIIE (4,90) (Z (J) , J= 1 , M) 

CA LL PROD (G, E, 1, M, 1 , G E , fi D , :1D , ID) 
C WHITE (4, 90) (GE(J) ,J=1,N) 



CALL ADD (XKKM1,GE,J;, 1 ,XKK,ND,MD) 



C 



CALL ?RCD (PHI,XKF,N,N, 1,XKKM1 ,yD,MD,LD) 
WRITE (0,8011) 

WRITE (4, 30 1 1 ) 
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3C 1 1 






C 

c 



L 

c 



38 13 



AT (/ ' 



') 



«RITE(J,30) (X^F (J) ,J = 1,R) 
mPITE {«, 30) (XKK (0) , J= 1 ,N) 
31 1 (K) =G (1, 1) 



GJ1 (?) =G (3, 1) 



?1 1 (K) =PKK«1 (1,1) 



?K1 1 (K) =EK? (1,1) 

Y3 (K) =XF'K ( 1) 

XH (?) =X?K (2) 

?TI’1E = IN1 (Ii:i3 (K) ) 

?/<TI:-1E = FlC AT (KTI'IE) 

1? (RKTl:i£. i^E.TIHE (?) ) GO'IO 88 13 
(?) =JULiiF 



EOB 1 (K) = YH (K) -LAI (X) 



EOB2 (X) =XH (X) -LONG (K) 



CONTINUE 



C WRITE(4,a8U) 

C UPITE (8, 83 12) 

8312 FC^UATI/' XKEUI ') 

L RITE (3, 9 0) (XKKU1 (J) ,J=1,N) 
'XRITE (4,90) (iXKUl (J) ,J=1,N) 

Y'HP (X) =XKKX1 (1) 

XHP (K) =XKKN 1 (2) 

HB=1300+ IT«100 
9685 DO 9634 1^1,4 

9036 IF (Ii:i (K) -£ 13. LE. 0) GC IC 99 37 
Et= JH+6 

IF (B5-II3E (K) . LI . 0) GC TC 36oo 
9997 IF (Ba-TIUE (K) . GI. 0) GO TC 9o99 
YKUO (X?K) = XKK (1) 

C 



XXM6 (KKK) =XKK (2) 

GO TO 9698 

9699 CC = bB- ( 1500 + 11 * 100) 
XFIi E (4,5o5o) CC, ba 
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(3, 5656) CC, 5K 

5cdo FC3V a: (5X, 'CC= ' , fig. 4, ' B3=',F10.4,* IQV* 

IF (CC-2 4. EC. ■)) GO 10 97C0 
IF (1-4. LI. G) GC rC 9o94 
I? (TiaE (K+ 1) -fiS. LI. 0) GC lO 9634 
9700 Ba=3E-24*1CC 

9694 IF (TI;1E (X+ 1) -B3.LE. 0) GC 10 9o34 
If {3B-Ii:^E (K) . EQ. 0) GC IC 9698 
DD'I = ED-lIi1E (K) 

WFITE (4, 38 U) 33,DDI, IT, K 
'.■;EITE{a,.3812)B3,DDl,IT, K 

38 12 F0?f-1AT(/,' FE= ',F10.4,' JDT= ',F10.4,' IT= 

1 ' K=',12) 

CALL PHIDEL ( DD T , N , L , A , 3 , PH I , DEL , D 1 , D 2 , ND , liD , I C) 
CALL ?FCD (PHI, XKK, i1 , i4 , 1 , XKKH6 , MD, JD, LD) 

XHITE (8,43 12) 

4312 fCP.HAT(y,' XKK.H6 

* B3 lIHE (X) •) 

XSITE (8,90) (XrvKF.o (J) , J = 1 ,3) ,38,113?. (K) 

WEITE (4, 43 12) 

WHITE (4,9 0) (XKK.'I6(J) ,J=1,N),2B,II;-1E(K) 

YKH6 (KFK) =XKr:.l6 (1) 

XK36 (ECXK) =XKK36 (2) 

Ob^a XEXK^KKK 
TI = I 

WHITE ( 4, 3-? 69) 

8969 FOx3dAT(/,' YK36 XEHo XKK RB 

1 I ') 

'..HITE (4, jO) Yr 36 (XKK) , XXH6 (FKK) ,XKKK,E2, LI 
KKK=KKX+ 1 
Bfl=B 3+6 
9634 CCNTIHDE 

WHITE (3 , 48 lx) 

W HITS (3, 90) (XFKH6 (J) , J = 1 , N ) , B B , 1 1 d E (K ) 
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IT" (a , '48 1^) 

■•;RIT5 (4,90) ( vvK.ib (J) ,J=1,9) (K) 

3 2 22 (:< + 1) - 150 0) * (.0 1) 

IS= {II ’'IE (K +2) - 1 500) (.01) 

C 

'/(BITE ( 4 , 9o9c) IT 
9696 FCBilAT (5X, ' IT= ',14) 

1=24*11 

C 

S=24*I3 

1 = TI.'1E (K + 1 ) - IT * 10 J - 1 50C fl 
S^riE.E (K + 2) -13* 100- 1500 + 3 
0 ’I = 3 - T 
0 3 = 6 

u «Hrin (3, 7865) o:< 

XRIIE (4,655) 

WRITE (3, £55) 

C DO 3022 1=1, N 

C WRITE (3, 90) (PBF.1 1 (1,0) ,0=1, ^i) 

030 2 2 ..RITE ( 4 , 90) (T.<KK1 (1,0) ,0=1,N) 

CALL ?HIDEL(DT,N,L,A,E,FHI,9EL,D1,D2,ND,0D,LD) 

C CALL T^ANS (D£1,N, 1,EEL1,NC,1C) 

C CALL PROD (CEL, DELI, N , 1 , )i, Q ,ND, .‘10,19) 

C C.ALL CCNSI (W (1 , 1) ,r, n,N,Q, ND,MD) 

C LDLLDLDLCLDLCLCI0LDLDLDIDLLDL9LCL91DLCIDLDIDLDLDLJL 
C 

C WE I T3 (4, 54 4) 

C KBITS (6,344) 

C DC 7 134 1=1, N 

C WRITE (4,90) (Q (1,0) ,0=1, N) 

C7134 J5?ITE (8, 90) (C ( 1 , 0 ) , 0= 1 , N ) 

C 

C *** HAS .3TE? OF FILTER REACHED POINT 25? 

IF ( (K-25) . N'E.O) GOTO 2222 
C 
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C *** .<G’J7I''IZ 10 PLACE ELLIPSE DAIA IN F^LE 

C 

THE 1 = - 50*ATAiN 1 ( 1 , 2) / (P^KEi 1 ( 1 , 1 ) - PKK .1 1 (2,2))) 

SU2X= (P?;K2 1 ( 1 , 1) ^-PKKHI (2, 2) ) /2. + PKK11 (1,2) 

1 /SIN (2.+THE1) 

SIG2 Y= (PKK.1 1(1,1) +Pfi'K:i 1(2,2)) /2.-PKK:1 1 (1,2) 

1 /SIN (2.=>^'IHE1) 

SX= ( (S IG2X) 5) *5 
S Y= ( (SIG2Y) **. 5) *5 
FI=3. 14159265/12 
CT=CCS (THE 1) 

S I = SIN (IHE 1) 

CC 1981 IELLIP=1,25 
XI=IELLIF 

XP (lELLI?) =SX^-C0S (?i*XI) =!'CT-SY*3Ii(?I*XI) *ST + XH{25) 
1? (lELLIP) =SX*C03 (PT + XI) ♦ST + SY^SIN (PT=i‘XI) *CT<-YH (25) 
1981 X3IIE (19, 1982) X? (lELLIP) ,Y? (lELLI?) 

1932 FCRHAT (2F1 5. 7) 

C *** END OF EILIPSE CALCOIATICM 

C 

r* 

2222 CCNTINUE 
C 

C wFIIE IKK STaFf 

'wFITE (8 , 2224) 

2224 fCr;XAT(//,' 1KN72(JJ), IKX48(JJ), rKH2'4(JJ) 

1 , IKli6(JJ) ') 

XEITE (8,2223) (TKd72 (JJ) ,T!<!>143 (JJ) , IK‘I24 (JJ) , IKX6 (JJ) 
1 ,JJ = 1,NN) 

WHITE (8,2225) 

2225 FOK.MAT (//, ' TKP72(;i1), TKr43(:iX), TXP24( -IM), 

1 TKP6(H:-1) ') 

WFITE (8, 2223) (TKP72 (MM) , '1/048 (/u1) ,TKP2 4 (NX) ,TK?t (NX) 
1 ,Ni'1=1,MN) 
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2 2^3 



YCR3AI (4110) 



C 

C '-.' 5 i I E { a , 7 7 7 7 ) 

i.’^IT2 (8,:>9 3) iTI'IZ (I) ,LAT (I) ,LO:iG (I) , YH (I) ,X[] (I) ,1 = 1 
1 ,iil) 

iVPITS (d, 39 j) {BTlaE{D ,ELAT (I) ,ELCi'iG{I) ,YKil6 (I) 

1 ,XK36 (I) , 1= 1 , NZ) 

393 FCP.MAT (U,5F 8.2) 

C 7<PITS (4,77 77) 

X HITE (4,393) (T 171 E { I ) , 1 A I ( I ) , LO .43 ( I ) , YH ( I ) , .XH (I) 

1 ,I=1,NZ) 

KaiTE (4, 393) (5TI3E (I) ,£IA‘T (1) ,aiC3G (I) , YKlifc (I) 

1 ,XX.'15 (I) ,I=1,NZ) 

c 

INTIi-!E = It1T (Tl:iE (K) ) 

PLTII’1E = FLCA1 (I4II3Z) 

If (ELI 171 E. 4E. ilYS (K) ) GGTC 399 

CALL T E L a ( E L A 1 , 3 L C N G , T I Z £ , H T I 3 E , X H , 1 a , L A I , L C M G ) 

3)9 CONTINUE 

C CRITE FOE1 AND EC32 TO FILE 'ECE DATA* 

K FITE (6, 32 9) 

329 f Oa:i AI (//, EX, ' JULIAN HR S03 1 

1 EO02’) 

XRrXE (6,323) (KJULHR (lECc) , SOB1 (lECB) ,ECS2 (IE3E) 

1 ,IECD=1,NZ) 

WRITE ( 1 6 , 3 28) (KJULHj (lECE) ,E0E1 (lECd) ,EC32 (lECc) 

1 ,IECD=1,NZ) 

328 rCF.dAT (I15,2F15.2) 

C 

DO 3^7 I=1,NN 

EY (1) =YK:16 (1) -BLAI (I) 

EX (I) =XKMo (I) -BLOMG (I) 

Ti:i£ (I) =ELCAT ( I) 
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j 27 Coi'iriM’JF 

;.TirE 3D J) (BTI^.E (I) , ELAT (I) , ELCNG (I) , Ei’ ( 
1 = 

'25112(8,38^0) (Gil (K) ,K= 1 , 10) 

OKITE (8, 5808) (TIc«2 (.3) 10) 

9898 FG5:iAT (!■ 10.4) 

’..’aiTE (3,4 10) 

4 10 FCSM AT ( ' 1 ' ) 

C n BITE (4, 7797) 

C 

C CALL PLGIP (iI8H,G1 1,NN, 1) 

C CALL BICIF (T i;'!£,G3 1 , LiN , 3) 

’.vBITE (8,410) 

C W H 1 1 E ( 4 , 7 7 7 7 ) 

C CALL PLOT!' (TI3E/ ^H,NN, 1) 

C CALL PLOIF (Ii;i2, iHF,AM,3) 

L L = ■ 1 + 2 

LONG (LI-1) =100 
LCNG(LL)= 160 
LA I (LL- 1 ) =0 
LAT (LL) =t0 
«HIT 2 (4 , 77 77) 

C CALL PLCII (LONG, LAI, LL, 1) 

C CALL PLCII (XH, )!H,MN,3) 

C CALL PLCTF (LONG, LAI, LL, 1) 

C CALL Pier? (XH, Y;i,:'iN,3) 

C 

c + BEITE VALUES I4TC FLCT FILE FOB DI3S!^i,A 

C MU.'IEEB CF VALUES 

'2 E IT E ( 10, 2 42 8) NN, 2 ( 1 , 1) ,5 (1,1) 

WPITE( 4,242a) N N , vi ( 1 , 1 ) , P ( 1 , 1 ) 

2 42 3 FORMAT (14, 3;< , ' W (1 , 1) = • , F3. 3, 3X, • R ( 1 , 1) = ' , F8 
C 

C CCL'JMN HEADINGS 



1) ,HX (I) 
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' 3 j 1 ' , c ' 



wFirZ (10,24^9) 

fcs:ia'l ('INT.Ea ' , 3 X, 'TI. 1 E' , 4 -(, 'Gil ' ,e;{, 
• YH? ' ,5i( , ’ LC^G • , 5X , ' LAT • ,7 a , ' Xn • ) 
C 



VALUES :0 EG PICIIHD 

WHITE (1J, 243C) {I,li;*!E <I) ,G1 1 (I) ,G31 (I) ,YH(I) /XnH (I) , 
*LCMG (I) , LAI (I) ,XH (I) , I = 1 , NN) 



2430 FCEHAI (I4,8FS.2) 



C 

WHITE (8, 10 16) A 
G « III IE (9, 10 16) K 

DO 65 0=1, N 

W HITE (8, 90) (G ( 0, I) , 1= 1 ,M) 

65 CONTINUE 
WHITE (3, 1o) 

16 FOnflAIC ',/,3X,'CCV. MAI. OF PBEDICIED ZSII.iAIE' ) 
DC 1 J 0= 1 , N 

WHITE (8, 90) (?KK:1 1 (1 , 0) , 1= 1 , N) 

13 CCN'TIMUE 
880 3TG? 

E2JD 



C W 3 W rt M I W 3 W .-1 W G 9 M M w M W A W M 'a M W A M W 3 W M W M W W M W M W M W M 3 w M W 3 V< M W 3 W d W W 3 W X 



SrijHOUTINE FillDSL (T , N, 3, A, E,PHI, DEL, D1 , 02, ND ,MC, ID) 
DIlENSIOIi A(12,12),8(12,12),?8I{12,12),DEL(12,12), 

♦ TEEM (12,12) , 

* COR ( 12, 12) ,C (12, 12) ,D1 (1 2, 12) ,B2 ( 12, 12) ,TE1L (12, 12) 
TEST = 1.E-7 

F=1 . 

DO 10 IR = 1,M 
DC 10 IC = 1 , N 
p[iI(IR,IC) = 0- 
?dI(IR,IF) = 1. 
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C(I?,IC) - A (I?,IC) 

I3JL (IR, IC.) = :/2.00*PuI {1H,IC) 

10 (I? , i>-) = (I? ,IC) 

5 0 DO 11 IE = 1 ,>I 

DO 11 IC = 1 

COR(IP,IC) = T/f*C(I5,IC) 

PHI(IR,IC) = PHI (IH,IC) +CCR (IH,IC) 

TEIL(IR,IC) = TZIl (Ib,IC) *-T/ ( (?+ 1. )*{F + 2. )) 

1 ■i'CCF (IE, IC) 

11 lEP:*! (IR, 1C) = TEEa (IP,IC) <-'1/ 1- )*CCK{IH,rC) 

DO 12 IR = 1,N 

DO 12 IC = 1 , N 
C(IR,IC) = 0. 

DO 12 K = 1,M 

12 C(IH,IC) = C (IS,IC) +A (IR,K) *COR (K,IC) 

F = F+1. 

DO 13 IE = 1,i1 
DO 13 IC = 1,N 

IF (APS (CCI; {IS,IC)).3T.IEST*ABS(?HI(IE,1C))) 

1 GO TO 50 

13 COWTii!10E 

CALL PRCC {TER1,B,a,L'l,a,DEL,HD,aD,LD) 

CALL PROD {TZIL,B,a,:i,;'l,D2,ND,flD,ID) 

DO 14 IG = 1,N 
DO 14 IC = 1 , a 

14 D1(IF,IC) = DEL (IH,IC) -D2 (IR,IC) 

PZIUoW 

5ND 

C Till3 SUdGOOTIME CO’-IPOTFS IHF CPTiaUM GAIM ’lAlBIX AND IH 

C CCVAEIANCZ 

C 

C 

C 

C 
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1 



K) 

ZKK { 12, 12) ,Q (12, 12) , :i (12, 12) , J (12, 12) 
1 ,2(12,12), 

1 ;-iIv12,12) (12,12) ,IF:^.?(12,12) ,I£2?1 (12,U) , 

1 (12,12) , 

6PKI { 12, 12) ,FHIT (12,12) , P:<K 11 (12,12) 

C G(K) = ? (K/K-1 ) (H*P (K/K- 1) *RT + 3) 

C PHI#? (K- 1/K- 1) ♦Frill + C 

IF (K. SQ. 1 ) GC TC ri9d9 

CALL T3AMS (PHI, N , N , FHIl, HD, HD) 

CALL P3CD (?K K, PHIT , N , N , N , IFHP, HD , HC, I D) 

CALL PHCC (PHI ,TEHP, H , N ,N , lEH? 1 , ri L' , HD , LD) 

CALL ADD (1 F:1 P 1 , 3 / N / / 2 K K ,1 1 , HD , .1 D ) 

WRITE (9, bS5) 

bbS FOHHAT(/' HAT3IX PKXMl •) 

DC 3022 1= 1 , N 

W3ITL (8,90) (PKK'II (I ,J) , J=1 , H) 

3022 WEI2E {^, (PKKM1 (I ,J) , 0=1 , N) 
dans CC N’T I HUE 

CALL T3ANG ( ri , H , N , HI , NU , H D) 

C t; 31 12(8,39) 

39 FOriHATC ii ') 

C DC 22 1=1, H 

C 2 2 WRITE (8, 90) (H (I, 0) , J= 1 , N) 

90 FGPHAT ( 1P7£ 1 1 . 4) 

C WRITE (3,3c) 

JO FoSHATC H'l ') 

C DC 23 1=1, N 

C 23 W RITE (8 , 90) (9T (I, J) , 0 = 1 , ;i) 

CALL PROD (PKKHI , HT,H, N, y ,IE:1F, HD, HD, ID) 

CALL PROD (R, TEH?,;-! , H, ;1, lEriPI ,ND,y:D, ID) 

C HRiTS(a,3d) 

38 FODHATC Ii ? HI •) 

C DC 24 1=1,. 1 
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C ^4 WFITE {3, 90) (TE:1P1 (I ,J) , J=1 , 

CALL .»Du (TE.'l? 1 [£3r 1 , L’D , ;1D) 

CALL 3 EC I? r , 0 . 0 JOCOO 1 , TE:1? 1 , 13:1 ?2 ,Kii3 , -iO) 
C IS i?2 ( 1 , 1) =is;-l?3 (2 , /DIZ 

C IE^1?2 (2, 1) =-'IEi12 3 (^, 1) / LSI 

C TF.;i?2 (1 ,2) =-IS.'123 (2, 1)/DSI 

C TE.1P2 (-t, 1) =-'LEi1?3 (1, 1)/D£I 

C WRITE (8,31) 

C DC 27 1=1,2 

C W3ITE (3, 90) (TE2P2 (I , J) , J= 1 , :i) 

3 1 rCR:1AT {' (K?H +5) - 1 ' ) 

IE {KES-2) 101,110,101 

1 10 WRITE (3,111) 

111 FCHiiAT (5[iKEF=2) 

101 CALL ?RCD (TE'1P,'iEM?2,N,2,M,G,i'iD,;iD,IC) 

C MOTE HERE Pi{K(I,J) = T (K/K) WHERE P(K/K) = 

C (I-G (K) *H) *P (E/K- 1) 

CALL ?FCD (G, h, N,.1, '1, I EE I , W D , 2£ , I C ) 



c 




WRITE (3, 30) 




c 


30 


FCR’IAI 1 ' 


GO ' ) 


c 




DC 25 1=1, N 




c 


25 


WRITE (3,90) 


(ISH?(1,J) ,J=1,W) 






DC 108 I = 


1 ,N 






DC 108 J = 


1,N 




1 Jd 


TE.'1P(I,J) = 


-te:i? (I, J) 


c 




WRITE (3, 37) 






37 


FCRMAI { • 


HI • ) 


c 




CO 45 1=1, W 




c 


4 5 


WRITE (3,90) 


( HI(I,J) ,J=1,N) 






CALL ADD (HI 


,TEi?iP, N, N ,TE2P, ND, .ID) 


c 




WRITE (8,33) 
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P(K/K-1) 





j3 


FCPNAT 


(' I 


-OH ') 




r' 

V- 




DC 35 


1=1, N 






C 


35 


W R 1 1 Z ( 




llEEP (I, J) ,J = 


1,N) 


C 




NCTE 


HERE 


PKK:-i1(I,J) = 


?(K/E-1) IlnEE E 



CALL 

-!5T!JPN 

LMD 



c 

c 

c 

c 





S'JiiPCUT 


IME 


ADD (A, 3, 


Im,K,C, ND,:-1D) 




DI/iENS ION 


A (ND,:iD) , 


o (MD, ,1J) ,C (ND,HD) 


L 


DO 2 


I = 






c 


DO 2 


J = 


1 




c 


2 C(I,J) 


= 0. 








DO 152 


I = 


1,N 






DC 152 


vJ = 


1,M 






15 2 C ( i , J ) 


= A {I , J) + E 


a, J) 




KZTU3N 









END 

C 

c 

r~ 

V- 

c 

SUHEC'JriNE iUD { A , E , N , .'I , C , ?.'D , D) 
DI‘1SitSI0I. A ( 'ID, -IE) , 3 (.\D, .ID) ,C (ND, .^C) 
DO 152 I = 1 
DO 152 J = 1,'1 

152 C (I, J) = A (I, J) - E (I , J) 

P.ETfJRN 

END 

C 

c 
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1 



151 



C 

C 

C 

C 



153 



C 



3;jri--vCuri:«F. prod ( a, l :jd, «c , lo) 

DI'IENSICN A (ED,MC) , 5 (?'. C, LD) ,C (ND,!^) 
DO 1 I = 1 , y 
DO 1 J = 1,L 
C (I, J) = 0. 

DO 151 I = 1,3 

DO 151 J = 1 ,L 

DO 151 K = 1,3 

C(I,J) =C(I,J) + A (1,K) *3 (K,J) 

HEIUSN 

EilD 



SOBROU nWE TRAMS (A,N,3,C, iND,3D) 
DIMENSION A (ND,MD) ,C (ED,:iD) 

BO 153 1 = 1,3 
DO 153 J = 1,3 
C(J,I) = A (I,J) 

RETURN 

END 



0 



c 



SUBROUTINE CONST ( ^ , A , N , 3 , C , N D , 3 D ) 
DINZNSIC N A ( ND, 3C} (MD , 30) 

IF(Q) 11,13,11 



10 DO roo I = 1,N 
DO 100 J ^ 1,3 
100 C (I, J) = 0.0 
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1 1 

12 

120 

13 

1 4 

14 0 

15 

150 

1 

2 

10 

1 1 

U 









IF (;-i 


. 


0) 13,U,13 


00 120 


1 


= 1,M 


DO 120 


J 


= 1,1 


0 (1,0) 


= 


A (1,0) 


EE rOFN 






IF (Q4l 


• 


0) 15,14,15 


DO 1 40 


T 

± 


= 1,N 


DO 1 40 


J 


= i,x 


C(I,J) 


- 


-A (1,0) 


HE'IUEO 






DO 150 


I 


= 1,0 


DO 150 


J 


= 1,0 


C (I, J) 


= 


C*A (1,0) 



3SIUHN 

£i1D 



S0 3KCOTIN2 ( N , £ ? , A , X , K E P , '1 ) 

DI.1EN3ION A ,X (0,i1) 

DO 1 1 = 1 , :i 

DO 1 0 = i,;-i 
X' ( I , J ) = 0 , 

DO 2 K = 1 , N 

= 1 . 

DO 34 L = 1,0 
K? = 0 
Z = 0. 

DO 1 ^ K = L , M 

I? (Z. 'JE. ABS (A (K , J) ) ) 00 10 12 

4 = AE5 (A (K, L) ) 

KP = K 
COOTIN'JE 
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13 

1 4 

15 

20 

30 

j 1 

32 

33 

35 

36 

34 

4 0 

4 1 

42 

43 

50 



1 ? (L.Or.KP) Go :c 23 
30 1 H J = L , 

Z = A (L, J) 

A (L, J) = A (KF, J) 

A (K?,J) = Z 

DO 15 J = 1 , -N 
Z = X{L,J) 

X (L, J) = X (KF,vl) 

X(KP,J) = Z 

IF (A33 (A iL,L) ) . LE, Z?) GC TO 50 

IF(L.3S.iN) GO TO 34 

L?1 = L+1 

DO 3 6 K = L F 1 , i1 

IF (A (K, L) . EC .0. ) GO TO 3o 

aAIIC = A {K, L) /A (L, L) 

DO 33 J ^ LP1,M 

A(K,J) = A (K, J) -BATIO^A (L, J) 

DO 35 J = 1,.\ 

X (X, J) = X (K,J) -HAIIO*X (L, J) 

Co ill I MiJE 
CO NT I MU E 
DC 43 1 = 1/N 
II = N+l-I 
DO 43 J = 1,N 
S = 0, 

IF(II.GE.K) GO IC 43 
IIPl = II + 1 

DO 42 K = II?1,M 
3 = 3 + A ( I 1 , K ) * X { K , J ) 

X(I1,J) = (X (II, J) -S)/A (11,11) 

2EIU0N 

KEP = 2 

SEIUFN 

E N D 
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c 

c 

c 

L 



c 

c 

c 

c 



c 

c 

c 



o U b K 0 D X _[ t i E *'i tv E A D { A , , Jv D , il D , 1 3 i A D ) 

DI'IEASIO.'; A (ND^'iC) ,IFE.AC (10) 

DO 10 I = 1,N 

10 HEAD (5, 20) (A (I, J) ,J = 1,M) 
zO fcp:iat (3f 10. 5) 
aEI'JFO 
E.ND 



SOBKCOTIHE MI.’hl’XE { A , N , « , tlD , .ID , I B I T £ ) 
DIOENSIC^ A (ND,''1C) (10) 

DO 10 I = 1 , N 

.’BITE (4,20) (I,J,A(I,J), J = I,'!!) 

10 WP.ITE(3,20) (I,J,A(I,J), J= 1,M) 

20 F0RHAaM2 (2X, ' (’ t 02, ’ , ' ,12, ♦) = ’,1?E10.3)) 

HETOEN 
END 



3 00 FO nil ME TE~ B (OLAi , BLCNG ETIHE , Al-i, YH, LAX 

1 , L 0 N V, ) 

DI lEMSION TIME (100) ,3TiaE (300) , ELAT (300) ,SLCNG (300) 
1 ,Yii(3C0) 

DIMELSION \U (300) , EDI (300) , EB2 (300) ,IBJUL (30 0) 

1 ,13 001(300), 

* JIIMS (300) , EOfal (300) , EC32 (jOO) 

PEAL + 4 L AT (300) , LONG (300) 
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r 



C 

o 

c 

c 



PJiEDTCIIuNS AND 



JULIAN Tl'tZ AND 
=*'■«'* 'JKSIIS DAIA'- 



?o l,i fCSIIIQN OF KALIiAN FIL'IE? 
BESI TRACK VALUES AS A FUNCTICN C 
NPIUS IdZ DATA TC THE FILE 



DO 26 1=1,50 
I3IIJE = INT (TINS (I) ) 

CALL JUTIiNE (ISIIZE,JULKE) 

I3JUL (I) =JUI. HE 
DO 10 J=1,50 
IBEST = INI (cTlJE (J) ) 

CALL jnilNE (IbZS'l, IBJ) 
IBJUL(J)=IEJ 

IF {IBJni (J) . NZ. ISJUl (I) ) GOTO 10 
EB1 (I) =YL (I) -bLAT (J) 

EB2 (I) =XLi (I) -BLCKG (J) 

C ZOB1 (I) =YH (1) -LAI (I) 

C 3032 {!) =XH (I) -LCN3 (I) 

E0B1 (1) =IAT (I) -FLAT (J) 

Z032 (I) =LCNG (I) -ELONG (J) 

JTi:i£ (I) =IEJ'JL (J) 

10 CONTINUE 

26 CONTINUE 

CGIIE (14,290) 

N 3IIE (o, 290) 



290 



1 



1 



rOE.lAT {//, 16X, • JIIJE' ,5X, ' EB1 • ,6X, ' E t 2 ' , 5 X , 
' ZOB 1 ' , 5X , ' ECB2 • ) 

DO 310 N=1,50 

IE (JTI.iZ (N) . EQ. 0) GOTO 310 

WRITE (6, 300) (JTI:‘:F(N) ,EB1 (N),EB2{N) ,2031 (N) 
,EC32 (N) ) 



.J RITE ( 1 4 , 300) (JTIJE (N) , E3 1 (N) ,EB2 (N) 



,203 1 (N) 
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1 ,E0d2C^j) 

300 FuOii .AT ( 1 5X, I 9, 4F9. 2) 

in CO'ITIOUZ 

C 

OSnFN 

cMD 

c 

c 

c 

SUBROailNii JUTIdE (IT 1 3E , J OlnK ) 

Q LCGNC 2C0IINE 

C CA.LCULATES JULIAN TIMS FPCd TEAS 1903 

IYF = AEA:w i;'jC = MCNTti (dAECH) , I£A = OAY, 

C Iti2 = HCUR OE lAY 

IYe=1982 
IdO=3 

I0A=I1I1E/ 100 
ina=ITI.1£-IDA* 100 

CALL YOdCSN (lYE , I dC , I C A , Id H , JO L HE ) 
r< 3 1 j 3 ;j 

c 

END 

SUBROUTINE N Q. iCS N ( Y E A R , dO , DA , H R , J UL d 3 ) 

C CALLFE BY SOERCUIINE JOIIdS 

C CALCULATES JULIAN CAY AND JULIAN HCUF 

INTEGEL IN ID (12) ,YE AR,DA, li? 

CAT A IN IC/0, 3 1,59, 9 0 ,120, 151, 181, ^12, 2 43, 273, 334 
1 ,334/ 

ID= (YEAF- 1900. ) «355.25-0.25 
IADD=0 

IF (dOD (YEAR, 4) .GT. G)GCTO b03 
IADD=1 

b03 JULDA = INID (dO) + DA + LADD 

IHR = 24. (ID + JOLCA - 1) + HF + 0.5 
IYEAF = YEAr.- 1 900 
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F. D 
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